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I  INTRODUCTION 


BACKGROUND 

This  study  of  the  water  quality  of  the  Oconee  River  Basin  was  conducted 
as  part  of  the  expanded  scope  Flood  Plain  Information  (XFPI)  study  done 
by  the  Savannah  District  Army  Corns  of  Engineers.  The  Savannah  District 
undertook  this  first  pilot  XFPI  study  In  which  qeographic  data  banks 
were  used  as  the  basis  for  simulating  watershed  hvdrology,  and  the  com¬ 
putation  of  expected  annual  flood  damages.  Environmental  considerations 
were  also  Included  in  the  original  pilot  study  objectives  but  mainly  for 
appraisal  of  wildlife  habitat  and  tradeoffs  between  the  desirability  of 
certain  land  uses.  The  primary  Intent  of  the  XFPI  analysis  was  to  analyze 
the  effects  of  alternative  futures  of  the  Oconee  River  Basin  development 
In  a  systematic  manner  such  that  realistic  comparisons  of  flood  hazards, 
flood  damages,  and  environmental  quality  could  be  made  between  existing 
and  alternative  future  watershed  develownent  patterns. 

After  the  viability  of  the  XFPI  geographic  data  bank  methodology  had  been 
successfully  demonstrated  for  flood  hazard  and  damage  computations,  the 
Savannah  District  requested  the  HEC  to  perform  an  Oconee  River  Basin 
water  quality  study  consistent  with  the  XFPI  objectives  and  methodology. 
Thus,  the  existing  and  alternative  future  water  quality  of  the  Oconee 
River,  within  the  study  area,  would  be  simulated  using  the  geographic 
data  bank  as  the  basis  for  land  use  inputs  to  existing  HEC  water  quality 
simulation  models.  The  Storaqe,  Treatment,  Overflow,  Runoff  Model  (STORM) 
[1]  would  be  used  for  determining  the  quantity  and  quality  of  land  surface 
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runoff  and  dry  weather  flow  and  the  Water  Quality  for  River/Reservoir 
Systems  Model  (WORRS)  [2]  would  be  used  to  simulate  water  quality  In  the 
river  network.  That  is,  the  land  surface  runoff  from  STORM  would  be  Input 
to  WQRRS  which  would  combine  all  Inflows  to  the  Oconee  River  and  simulate 
the  resultant  river  water  quality. 

Some  historical  data  about  the  water  quality  of  the  Oconee  River  Basin 
were  said  to  be  available  but  the  extent  and  appropriateness  of  those  data 
for  these  modeling  purposes  were  not  known.  The  HEC  undertook  the  water 
quality  studies  expecting  that  the  historical  data  would  be  satisfactory 
for  model  calibration  for  existing  conditions.  If  some  aspects  of  these 
data  were  not  sufficient,  then  general  experience  from  other  water  quality 
studies  would  be  used  to  ascertain  acceptable  performance  of  the  simulation 
models. 

The  Savannah  District  also  requested  the  HEC  to  study  the  land  surface 
erosion  simulation  possibilities  using  the  geographic  data  bank.  In 
particular.  It  was  desirable  to  erode  and  transport  sediment  on  a  grid 
cell  basis  as  defined  by  the  topography  stored  In  the  geographic  data  bank. 
A  new  computer  program  would  be  developed  to  Implement  this  proposed 
methodology. 


SCOPE  AND  OBJECTIVES 

The  objective  of  the  HEC  study  was  to  Investigate  the  applicability  of  the 
HEC  water  quality  simulation  models,  STORM  and  WQRRS,  for  usage  in  XFPI 
studies.  This  was  to  be  accomplished  through  an  evaluation  of  the  water 
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quality  Impacts  of  existing  and  alternative  future  land  use  development 
patterns  In  the  Oconee  River  Basin.  The  new  HEC  grid  cell  sediment  transport 
model  would  be  evaluated  In  a  similar  manner.  The  methodologies  for  analysis 
of  water  quality  and  sediment  transport  were  to  be  consistent  with  the 
philosophy  (l.e.,  use  of  geographic  data  banks)  of  the  ongoing  XFPI  pilot 
study. 

The  land  use  data  required  by  the  STORM  model  were  to  be  obtained  from  the 
geographic  data  bank  through  the  Hydrologic  Parameters  (HYDPAR)  utility 
program.  STORM  would  also  access  the  land  use  data  for  computation  of 
sanitary  sewage  flows.  Certain  changes  were  reouired  In  STORM  to  utilize 
the  grid  cell  data  bank.  The  HEC  would  calibrate  the  models  on  existing 
data  and  use  the  calibrated  models  to  simulate  alternative  future  watershed 
developments. 

STUDY  TEAM 

This  study  was  to  be  carried  out  entirely  by  the  HEC  with  minimal  direct 
Involvement  of  Savannah  District  personnel.  The  district  was  to  provide 
general  guidance  about  the  XFPI  study,  the  objectives  of  this  study,  and 
supply  existing  water  quality  and  sediment  data  to  the  HEC.  Because  this 
study  was  a  special  investigation  of  the  applicability  of  new  techniques 
in  support  of  the  XFPI  study,  the  district  did  not  feel  it  was  necessary 
to  have  this  water  quality  modeling  expertise  developed  within  their  staff. 
The  district  staff  could  be  trained  in  the  use  of  this  methodology  at  a 
later  date  if  warranted  by  their  needs  and  the  results  of  this  study. 
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The  HEC  conducted  the  study  as  a  team  effort  with  Mr.  Jess  Abbott  being 
responsible  for  the  application  of  the  STORM  model.  Dr.  Michael  Gee  being 
responsible  for  the  grid  cell  sediment  transport  Investigations,  and  Mr. 

R.  G.  Willey  being  responsible  for  the  application  of  the  WORRS  model. 
Messrs.  Darryl  Davis  and  Pat  Webb  provided  guidance  on  the  XFPI  methodology 
and  the  utilization  of  the  geographic  data  bank.  Mrs.  Marilyn  Hurst  and 
Mr.  Paul  Ely  performed  most  of  the  detailed  tasks  involved  to  complete  the 
project.  The  drafting  was  done  by  Mr.  Roger  Nutter. 
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II  SUMMARY  AND  CONCLUSIONS 


SUMMARY 


The  study  objectives  were  carried  out  as  proposed.  The  development, 
modification,  and  implementation  of  the  mathematical  models  provided  the 
means  to  simulate  exlstlnq  and  assess  the  future  water  quality  and  sediment 
transport  characteristics  of  the  Oconee  River  Basin  study  area.  The  STORM 
and  WQRRS  water  quality  modeling  methodology  was  successfully  Implemented; 
however,  the  lack  of  adequate  data  to  calibrate  these  models  made  their 
application  difficult  and  the  results  unsubstantiated  except  for  general 
comparisons  with  experience  from  similar  studies. 


The  grid  cell  land  surface  erosion  and  sediment  transport  methodology  was 
not  entirely  successful  because  of  problems  In  using  the  topography  file 
of  the  geographic  data  bank.  The  erosion/sediment  transport  methodology 
was  based  on  continuously  downhill  sloping  grid  cells  to  the  stream 
collection  network.  Upon  application  of  the  method  using  the  Oconee 
topography  file,  it  was  learned  that  the  grid  cells  did  not  slope  con¬ 
tinuously  downhill  to  the  collection  channels.  A  concentrated  effort  was 
made  to  incorporate  the  required  slope  continuity  into  the  geographic  data 
bank.  This  problem  could  not  be  resolved  within  the  scope  of  this 
investigational  effort.  If  further  development  of  this  methodology  Is 
desired,  the  erosion/sediment  transport  methodology  must  be  modified  to 
accept  the  existing  data  bank  topography  or  the  topographic  data  must 
be  edited  to  conform  to  the  slope-continuity  assumptions  made  In  the 
mathematical  model. 


The  STORM  and  WQRRS  Interfaces  with  the  data  bank  were  developed 
simultaneously  with  assembly  of  historical  data  required  to  operate  and 
calibrate  the  models.  It  soon  became  evident  that  very  few  historical 
water  quality  measurements  had  been  made  in  the  Oconee  River  and  its 
tributaries.  Thus,  the  calibration  of  the  models  would  have  to  be  based  on 
theory  and  experience.  Some  water  auality  measurements  had  been  made  and 
these  data  were  used  as  much  as  possible.  No  data  were  available  for  the 
important  storm  runoff  periods  durlnq  which  land  surface  pollutant  washoff 
occurs.  Concentrations  of  pesticides,  heavy  metals  and  other  parameters  were 
not  specifically  mentioned  and  not  evaluated  in  this  study.  Such  parameters 
cannot  be  simulated  presently  by  STORM  and  WQRRS. 

The  decision  was  made  to  continue  with  the  STORM-WORRS  modeling  effort  with¬ 
out  adequate  data.  This  was  done  mainly  so  that  the  general  methodology 
could  be  developed.  Had  it  not  been  for  the  desire  to  develop-and  demonstrate 
the  STORM-WQRRS  methodology  for  XFPI  studies,  a  more  simple  water  quality 
study  would  have  been  recommended  for  the  Oconee  study.  The  recoirmended 
water  quality  study  methodology  would  have  been  commensurate  with  the 
detail  of  existing  data  and  the  degree  of  detail  warranted  by  the  study 
objectives. 

Assumptions  were  made  about  the  basic  water  quality  Inputs  from  the  Upper 
Oconee  River  (outside  the  study  area)  and  about  the  land  surface  runoff 
within  the  study  area.  Because  there  were  no  data  on  the  river  water  quality 
during  storm  events,  the  calibration  of  reasonable  values  of  land  surface 
runoff  and  Incoming  river  water  quality  required  much  more  time  than 
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anticipated.  The  land  surface  runoff  simulation  was  reviewed  In  some  detail 
for  Reaches  1  and  2,  Figure  III-3.  Upon  achieving  reasonable  results  for  these 
reaches.  It  was  assumed  that  subsequent  reaches  simulated  in  a  similar  manner 
would  also  have  reasonable  results.  Thus,  the  other  basins  were  simulated 
and  it  was  not  until  all  of  the  basins  were  aggregated  in  the  WQRRS  receiving 
water  model  that  it  was  apoarent  some  basin  results  were  unreasonable.  The 
complexity  of  the  WORRS  model  did  not  facilitate  the  timely  appraisal  of  these 
potential  problems.  At  that  late  date  in  the  study,  some  of  the  basin  land 
surface  runoff  simulations  had  to  be  rerun  and  the  river  system  was  simulated 
again.  The  new  results  were  acceptable.  The  land  surface  runoff  from  the 
basins  should  have  been  reviewed  more  thoroughly  during  the  initial  simulations; 
however,  such  review  was  not  emphasized  because  of  the  data  limitations. 

The  existing  land  use  condition  was  simulated  using  non-point  source  land 
surface  runoff,  point  sources  within  the  tributary  subbasins,  and  point  source 
Inputs  from  the  two  main  sewage  treatment  plants  in  the  study  area.  The 
treatment  plant  loads  appear  to  cause  the  most  significant  Impact  on  the 
North  Oconee  River  for  existing  conditions.  The  pollutant  loadings  from 
the  sewage  treatment  plants  on  the  North  Oconee  contribute  approximately 
80-93%  (depending  on  the  parameter)  of  the  total  loads  (point  and  non-point 
sources).  Therefore  further  reduction  in  the  loads  from  sewaqe  effluent 
would  appear  to  have  the  greatest  effect  in  upgrading  the  water  quality  for 
existing  conditions.  There  does  not.  appear  to  be  a  significant  water 
quality  problem  under  existing  conditions  because,  with  the  possible  ex¬ 
ception  of  conforms,  established  water  quality  standards  are  not  exceeded. 


The  basins  having  the  most  significant  Impact  are  1A,  2,  and  3  on  the  North 
Oconee  and  6A,  6B,  and  16  on  the  Middle  Oconee,  Figures  III-l  and  II 1-2. 

The  sewage  treatment  plant  effluents  on  the  Middle  Oconee  contribute  48-80% 
of  the  total  loads  reaching  the  Middle  Oconee  above  the  confluence.  The 
concentrations  do  not  exceed  established  water  quality  standards  under  the 
existing  development  with  the  exception  of  coliforms.  The  most  irinedlate 
Improvements  could  be  made  by  reducing  the  loads  from  treatment  plants. 

The  North  Oconee  River  watershed  contributes  a  significantly  larger  pollutant 
load  to  the  Main  Oconee  River  than  does  the  Middle  Oconee  River.  The  North 
Oconee  watershed  has  significantly  larger  loadings  from  the  sewage  treatment 
plants.  Specifically  the  North  Oconee  plant  contributes  about  70%  of  the 
total  load  (point  and  non-point  sources)  to  the  Main  Oconee  River.  In 
addition,  effects  of  the  loadings  from  the  treatment  Diants  are  much  more 
pronounced  In  the  North  Oconee  than  the  Middle  because  the  natural  flows 
from  the  upstream  watershed  on  the  North  Oconee  provide  much  less  dilution. 

The  future  water  quality  for  the  Oconee  River  study  area  was  simulated  for 
the  Alternative  C,  Table  III-l,  land  use  development  plan.  In  general,  the 
sources  of  water  degradation  are  the  same  as  those  defined  for  existing 
conditions.  On  the  North  Oconee  the  contributions  due  to  treatment  plant 
loadings  ranged  from  81-88%  of  the  total  loadings  and  on  the  Middle  Oconee 
the  plant  loadings  contribute  48-88%  of  the  totals. 

i 

j 

While  the  percentage  contribution  from  sewage  treatment  loads  appears  to  ] 

remain  constant,  the  total  loads  have  increased  somewhat  from  existing  to 
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the  alternative  future  C.  The  increases  are  not  major  due  to  the  relatively 
small  percentage  change  in  land  use  change  in  the  total  study  area.  The 
major  Impacts  were  shown  to  be  In  subbasins  1A,  2,  and  3  on  the  North  Oconee 
and  6A,  6B,  and  16  on  the  Middle  Oconee  because  these  subbasins  experienced 
the  greatest  degree  of  urbanization.  These  subbasins  and  the  loadings  from 
the  2  main  treatment  plants  tend  to  create  "shock  loadings"  In  the  reaches 
immediately  downstream  of  the  effluent  outfalls  or  tributary  Inflows. 


9 


CONCLUSIONS 


Conclusions  with  respect  to  both  the  technical  feasibility  and  the 
suitability  of  the  methodology  for  water  quality  and  land  surface  erosion/ 
sediment  transport  studies  in  support  of  the  XFPI  program  will  be  made  in 
this  section.  The  impact  of  the  future  land  use  plans  on  Oconee  River 
water  quality  were  discussed  in  the  preceding  Summary  Section.  The  study 
objectives,  budget,  and  availability  of  data  to  support  the  proposed  study 
methodology  are  important  factors  in  determining  the  most  appropriate 
technology  for  a  project.  The  STORM  and  WQRRS  water  quality  simulation 
models  have  been  shown  to  be  technically  feasible  for  water  oualitv  studies. 
The  appropriateness  of  their  use  for  XFPI  studies  is  a  much  more  important 
question  which  can  only  be  determined  by  the  water  quality  objectives  of 
each  study  application.  The  STORM  and  WORRS  computer  programs  and  attendant 
study  methodoloqy  are  quite  appropriate  for  detailed  water  quality  studies. 
That  is,  these  models  provide  a  good  simulation  of  the  ohysical  water 
quality  svstem,  both  for  land  surface  runoff  and  receiving  waters.  The 
STORM  model  provides  a  relatively  simple  simulation  of  land  surface  runoff. 
The  WQRRS  model  performs  a  rather  complex  simulation  of  receiving  water 
quality  and  requires  much  more  comprehensive  input  data  than  does  STORM. 

The  general  objective  of  XFPI  studies  Is  the  analysis  of  the  hvdroloqic, 
economic,  and  environmental  impact  of  future  land  use  development  patterns. 
To  accomplish  this,  the  existing  system  must  be  reoresented  satisfactorily 
in  simulation  models.  There  must  also  be  a  consistent,  logical  means  to 
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generate  and  compare  alternative  futures.  Because  there  are  so  many  unknowns 
with  respect  to  specific  location  and  type  of  future  land  use  patterns,  river 
regulation,  and  waste  water  management  facilities,  the  analysis  of  futures 
can  be  less  detailed  than  known  conditions.  Methods  should  capture  the 
esseoce  of  the  future  conditions  without  being  overly  complex  about  the 
specific  types  and  locations  of  the  development. 

In  accordance  with  the  above  objectives,  the  STORM  model  provides  both  the 
type  of  Information  and  the  level  of  technical  detail  which  are  appropriate 
for  XFPI  studies.  The  basic  land  use  parameters  of  the  STORM  model  are 
readily  derived  from  the  geographic  data  bank.  Other  Input  to  the  STORM 
model  can  be  easily  obtained  or  estimated  from  previous  experience. 

For  application  of  the  STORM  model,  continuous  rainfall  and  runoff  data  are 
recommended  for  a  period  of  several  years  to  calibrate  the  hydrologic  parameters 
of  the  model.  Pollutant  loadings  at  the  subbasin  outlets  should  be  measured 
throughout  several  major  storm  events  during  the  multiyear  hydrologic  calibration 
period.  If  these  data  are  not  available,  the  STORM  model  should  not  be  used 
unless  acceptable  results  can  be  obtained  from  use  of  coefficients  derived  from 
similar  studies.  The  availability  of  data  should  be  determined  early  in  the 
study  so  that  data  collection  efforts  cannbe  arranged  as  necessary. 

The  water  quality  simulation  capability  of  the  WQRRS  model  seems  to  be  more 
comprehensive  than  required  for  the  general  XFPI  study.  The  many  data 
requirements  limit  the  utility  of  this  model.  The  WQRRS  model  would  be 
required  if  a  more  comprehensive  understanding  of  the  water  quality  condition 
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Is  desired.  This  might  be  the  case  for  reservoir  regulation  studies  or 
major  river/reservoir  studies  with  specific  water  quality  objectives.  For 
the  XFPI  level  of  complexity,  a  more  simple  receiving  water  analysis  seems 
appropriate.  An  analysis  conmensurate  with  the  complexity  of  the  STORM 
model  Is  desirable.  The  HEC  Is  presently  developing  that  type  of  simplified 
receiving  water  model. 

If  a  detailed  water  quality  study  had  been  required  for  the  Oconee  XFPI  study, 
then  data  collection  efforts  should  have  been  started  as  soon  as  It  was 
determined  that  the  historical  data  were  Inadequate.  The  collection  and 
analysis  of  field  data  would  have  required  a  considerably  larger  study  budqet, 
on  the  order  of  8  to  10  times  the  Initial  water  Quality  study  budget  of 
$20,000. 

For  application  of  the  WORRS  model,  climatic,  pollutant  Dolnt-source  loading 
data,  and/or  results  from  STORM  must  be  known  for  the  entire  calibration 
period.  Most  Importantly,  in-stream  water  quality  measurements  must  be 
available  for  several  major  parameters  both  during  storm  runoff  periods 
(preferably  the  same  periods  as  used  for  STORM)  and  low  flow  or  other  critical 
water  quality  periods.  These  data  should  be  available  for  at  least  one 
location  in  the  river  system  depending  upon  how  much  variation  there  Is  In 
the  land  use  and  stream  regimes  in  the  river  network.  If  only  one  location 
were  available.  It  should  be  at  the  downstream  boundary  of  the  river  network 
so  that  the  Integrated  effects  of  the  land  surface  runoff  and  In-stream 
quality  changes  are  measured. 
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Ill  OCONEE  RIVER  SYSTEM  DATA  AVAILABILITY 
GENERAL 

The  Oconee  River  begins  in  the  Georgia  counties  of  Barrow  and  Jackson  north 
of  Athens  and  flows  south  through  the  middle  of  Georgia.  After  It  joins  the 
Ocmulgee  River  near  Hazlehurst,  it  becomes  the  Altamaha  River  which  flows 
southeast  to  the  Atlantic  Ocean.  A  location  map  is  shown  in  Figure  I I I -1 . 

The  study  boundaries  for  this  project  include  the  Oconee  River  drainage  be¬ 
tween  the  Currey  Creek  dam  site  on  the  North  Oconee  and  the  State  Highway 
Bridge  33  on  the  Middle  Oconee  down  to  a  location  8  miles  below  Barnett 
Shoals  Dam  on  the  Main  Oconee  (i.e.,  inflow  to  Wallace  Reservoir).  The 
Oconee  River  study  area  is  shown  in  Figure  1 1 1-2  and  schematics  of  the  study 
area  are  shown  in  Figures  1 1 1 -3  and  III -4 .  The  schematics  include  the 
location  of  all  major  tributaries,  sewage  treatment  plant  effluents,  and 
the  Athens  water  supply  intake. 

The  historical  period  to  be  used  for  analysis  was  selected  using  the 
following  criteria: 

(1 )  A  low  flow  period. 

(2)  A  period  with  several  significant  rain  events. 

(3)  A  relatively  recent  period  (i.e.,  existing  conditions). 

(4)  A  short  duration  (i.e.,  one  month). 

Water  quality  sampling  points  were  located  on  the  Middle  Oconee  River  in 
1970  and  on  the  North  Oconee  in  1974.  During  1970-1975,  the  dryest  one 
month  period  having  3  to  4  significant  rain  events  was  October  1970.  This 
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period  has  as  much  available  water  quality  data  as  any  period  in  1970  to 
1975  and  was  therefore  selected  for  analysis. 

The  various  aspects  of  specific  data  availability  will  be  discussed  in  the 
remainder  of  this  chapter. 


NORTH  OCONEE 
BOUNDARY  INFLOW 


Subbasin  20  (34.8 


Subbosln  17-Shonkles  Cr.  (33.5) 


Subbasin  18-  Crooked  Cr.  (31.8) 


Subbasin  19(31.0) 


Subbasin  lB-Cubb  Cr.  (28.7) 


Subbasin  1A  (26.0 


Subbasin  2-Sandy  Cr.  (24.3 


Athens  Water  Intake  (23.7 


Subbasin  3  (23.2 


Subbasin  4  -Trail  Cr.  (21.7) 


Subbasin  12B  (20.3 


North  Oconee  STP  (19.9) 


Sub  basin  12C  (19.0) 


Fig  HI- 3.  Schernatic  of  North  Oconee  River 
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MIDDLE  OCONEE 
BOUNDARY  INFLOW 
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Subbasin  5C  (31.2) 


Subbosin  5D  — Bear  Cr.  (29.5) 

Subbosln  5A—  Turkey  Cr.(25.9) 

Subbosin  5B(24.6) 

Middle  Oconee  USGS  Gage  2-2175(22.3] 
Subbosin  6B(22.Q) 

Subbosin  6A  ( 19.9) 

Middle  Oconee  STP(17.3) 


ii 


Middle  and  Main  Oconee  River 


Subbasin  7  —  Barber  Cr.  (17.0 
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Subbasin  15  (9.8) 


Fig.  m-4.  Schematic  of  Middle  and  Main  Oconee  River 


METEOROLOGY 


The  weather  data  for  the  analysis  were  obtained  from  the  National  Weather 
Service  using  the  Athens  Municipal  Airport  Weather  Station.  A  magnetic 
tape  of  hourly  rainfall  at  five  stations  for  the  STORM  model  input  and 
another  tape  of  dry  and  wet  bulb  air  temperature,  barometric  pressure, 
wind  speed  and  cloud  cover  for  the  WORRS  model  input  were  obtained  from 
the  Ashville,  North  Carolina  office  of  the  U.S.  Weather  Service. 

These  input  data  were  the  easiest  to  obtain  of  all  the  required  input 
data  for  either  STORM  or  WQRRS.  See  Chapter  IV  for  discussion  of  models. 
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LAND  USE 


Land  use  Is  one  of  the  most  Important  input  variables  for  STORM.  It  is 
especially  important  in  this  study  since  one  of  the  main  objectives  is  to 
assess  the  impact  of  future  development  (as  characterized  by  land  use)  on 
the  water  quality  of  the  Oconee  River.  Land  use  for  each  STORM  watershed 
was  taken  directly  from  output  from  HYDPAR  (an  HEC-developed  utility  program 
to  calculate  hydrologic  parameters  from  a  grid  cell  data  ban!;).  The  specific 
land  use  categories  that  were  used  in  this  study  are  as  follows: 


Code  No. 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Designation 

Developed  Open  Space  (lawns,  parks,  golf  courses, 
cemeteries  and  rights-of-way) 

Low  Density  Residential 

Medium  Density  Residential 

High  Density  Residential 

Agriculture  (cultivated  land,  row  crops,  small  grain) 
Industrial 

Commercial  (strip  and  isolated  commercial) 

Pasture 

High  Density  Commercial  (downtown  areas  and  shopping  centers) 

Institutional 

Natural 


Some  alterations  were  made  to  a  few  of  the  land  use  categories.  Low 
Density  Residential  and  Medium  Density  Residential  were  combined  into  a 
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single  category.  Hardwoods,  pines,  and  wetlands  were  combined  Into  a  single 
category  (Natural).  Roads,  land  fills  and  water  bodies  were  not  simulated, 
however,  the  areas  of  each  of  these  categories  were  subtracted  from  each 
watershed  area.  Table  III-l  shows  the  land  use  for  each  STORM  watershed 
for  both  existing  and  one  alternative  future  (1990C). 


TABLE  m-1 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

NO. 

SUB -BAS  IN 
NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

* 

k> 

mm 

1 

20 

1870 

Low  Residential 

1.1 

1.1 

High  Residential 

0.0 

0.0 

Commercial 

0.0 

0.0 

Industrial 

0.0 

0.0 

Pasture 

25.5 

25.5 

Natural 

73.4 

73.4 

1 

17 

3700 

Low  Residential 

2.3 

2.3 

High  Residential 

0.0 

0.0 

Conmerclal 

0.0 

0.0 

Industrial 

0.0 

0.0 

Pasture 

38.4 

38.4 

Natural 

59.3 

59.3 

I 

18 

5620 

Low  Residential 

1.8 

1.8 

High  Residential 

0.0 

0.0 

Commercial 

0.0 

0.0 

Industrial 

0.0 

0.0 

Agricultural 

0.8 

0.8 

Pasture 

52.4 

52.4 

Natural 

45.0 

45.0 

1 

19 

2180 

Low  Residential 

1.1 

1.1 

Medium  Residential 

2.4 

2.4 

Commercial 

0.0 

0.0 

Industrial 

0.0 

0.0 

Pasture 

33.9 

33.9 

Natural 

62.6 

62.6 

1 

IB 

3741 

Low-Medium  Residential 

2.3 

6.2 

High  Residential 

0.4 

0.4 

Conmerclal 

0.3 

0.6 

Industrial 

0.0 

0.0 

Agricultural 

31.3 

28.2 

Pasture 

7.5 

7.4 

Natural 

57.8 

54.8 

Open 

0.4 

0.6 

NOTE:  Water  Area  Is  Included  In  Natural 
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TABLE  III-l 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

l 

1 

1A 

2928 

Low-Medium  Residential 

11.3 

13.0 

Hlqh  Residential 

0.6 

1.8 

Commercial 

0.2 

0.6 

Industrial 

4.9 

11.2 

Agricultural 

9.7 

6.3 

Pasture 

11.1 

7.9 

Institutional 

0.2 

2.8 

Natural 

60.0 

51.3 

Open 

2.0 

5.0 

1 

2 

41254 

Low-Medium  Residential 

2.1 

3.1 

High  Residential 

0.1 

0.2 

Commercial 

0.1 

0.2 

Industrial 

0.0 

0.0 

Pasture 

41.8 

41.7 

Institutional 

0.1 

0.1 

Natural 

50.9 

50.0 

Open 

0.3 

0.4 

2 

3 

2272 

Low-Medium  Residential 

30.9 

33.1 

High  Residential 

5.9 

4.4 

Commercial 

4.3 

4.6 

Industrial 

14.4 

15.1 

Agricultural 

2.7 

2.1 

Pasture 

3.1 

2.3 

Hlqh  Commercial 

1.0 

1  .0 

Institutional 

1.6 

1.5 

Natural 

28.8 

26.6 

Open 

7.3 

6.7 

Roads 

0.0 

2.6 
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TABLE  III-l 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

W 

to 

1990C 

% 

2 

4 

7915 

Low-Medium  Residential 

20.0 

21.5 

High  Residential 

2.6 

2.6 

Corrmercial 

1.8 

2.0 

Industrial 

1.5 

15.3 

Agricultural 

24.0 

17.9 

Pasture 

6.3 

4.1 

Hiqh  Commercial 

0.0 

0.5 

Institutional 

0.7 

1.0 

Natural 

43.1 

31.8 

Open 

0.0 

0.6 

Roads 

0.0 

2.9 

2 

12B 

1756 

Low-Medium  Residential 

18.6 

18.6 

Hiqh  Residential 

4.q 

6.0 

Commercial 

4.2 

3.8 

Industrial 

o.n 

0.0 

Agricultural 

12.3 

9.7 

Pasture 

0.6 

0.3 

High  Commercial 

3.5 

3.5 

Institutional 

34.6 

34.1 

Natural 

1Q.9 

18.7 

Open 

1.4 

2.7 

Roads 

9.0 

2.7 
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TABLE  III-l 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

? 

1990C 

l 

2 

12C 

2702 

Low-Medium  Residential 

17.7 

19.7 

High  Residential 

3.2 

4.6 

Commercial 

3.8 

3.7 

Industrial 

2.4 

15.2 

Agricultural 

8.2 

5.3 

Pasture 

10.0 

5.5 

High  Commercial 

1.1 

0.9 

Institutional 

3.2 

2.8 

Natural 

39.0 

31.5 

Open 

11.4 

7.2 

Roads 

0.0 

3.5 

2 

12A 

3191 

Low-Medium  Residential 

7.9 

9.0 

High  Residential 

0.5 

2.8 

Commercial 

1.3 

1.3 

Industrial 

0.3 

1.9 

Agricultural 

11.8 

9.0 

Pasture 

5.1 

3.7 

High  Commercial 

0.0 

0.7 

Institutional 

3.4 

5.8 

Natural 

65.9 

60.5 

Open 

3.8 

4.8 

Roads 

0.0 

0.4 

3 

5C 

3456 

Low  Residential 

1.6 

1.5 

High  Residential 

0.4 

0.4 

Commercial 

0.5 

0.7 

Industrial 

0.0 

0.0 

Agricultural 

16.6 

16.5 

Pasture 

4.0 

4  .0 

Natural 

75.9 

75.9 

Open 

1.0 

-  -  -  _ 

1.0 

TABLE  III-l 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

SUB-BASIN 

AREA 

LAND  USE 

EXIST 

NO. 

NO. 

(ACRES) 

* 

3 

50 

13020 

Low  Residential 

2.3 

Medium  Residential 

1.2 

Commercial 

0.3 

Industrial 

0.0 

Agricultural 

0.9 

Pasture 

50.3 

Institutional 

0.2 

Natural 

43.9 

Open 

0.9 

3 

5A 

5723 

Low-Medium  Residential 

18.3 

High  Residential 

1.7 

Commercial 

0.6 

Industrial 

1.1 

Agricultural 

9.5 

Pasture 

3.9 

High  Commercial 

0.2 

Institutional 

0.6 

Natural 

62.8 

Open 

1.3 

3 

5B 

3385 

Low-Medium  Residential 

8.1 

Hlqh  Residential 

0.9 

Commercial 

1.1 

Industrial 

0.0 

Agricultural 

17.5 

Pasture 

4.7 

High  Commercial 

0.2 

Natural 

67.4 

Open 

0.1 

3 

6B 

2658 

Low-Medium  Residential 

23.9 

High  Residential 

5.4 

Commercial 

7.2 

Industrial 

0.0 

Agricultural 

0.9 

Pasture 

3.0 

High  Commercial 

1.7 

Institutional 

2.7 

Natural 

53.4 

Open 

1.8 

TABLE  III-l 


LAND  Uc-E  BY  WATERSHED  AND  REACH 


It- 

I? 


• 


f; 
*  . 


* 

5 

-4 

»  v 

fc 

<: 

♦ 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

* 

* 

1990C 

f 

3 

6A 

3443 

Low-Medium  Residential 

47.6 

49.5 

High  Residential 

2.4 

3.2 

Commercial 

5.9 

5.9 

Industrial 

0.0 

0.0 

Agricultural 

1.9 

0.9 

Pasture 

4.9 

3.7 

High  Commercial 

1.7 

1.7 

Institutional 

3.0 

2.8 

Natural 

32.0 

29.5 

Open 

0.6 

1.2 

Roads 

0.0 

1.6 

3 

7 

27410 

Low  Residential 

3.0 

2.8 

Medium  Residential 

0.8 

1.5 

Commercial 

0.1 

0.2 

Industrial 

0.0 

0.2 

Agricul tural 

16.0 

15.4 

Pasture 

37.1 

37.3 

Institutional 

0.1 

0.2 

Natural 

42.8 

42.3 

Open 

0.1 

0.0 

4 

8 

10260 

Low-Medium  Residential 

13.1 

15.2 

High  Residential 

0.8 

1  .4 

Commercial 

1.1 

1.2 

Industrial 

0.0 

0.0 

Agricultural 

16.6 

14.7 

Pasture 

24.8 

25.5 

High  Conmercial 

0.0 

0.1 

Institutional 

0.2 

0.5 

Natural 

43.0 

40.5 

Open 

0.4 

0.7 

Roads 

0.0 

0.2 

TABLE  1 1 1-1 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

NO. 

SUB-BASIN 

NO. 

AREA 

(ACRES) 

LAND  USE 

EXIST 

to 

H 

4 

9 

1290 

Low  Residential 

0.8 

0.8 

Medium  Residential 

8.3 

7.8 

Commercial 

1.2 

1.5 

Industrial 

0.0 

0.0 

High  Residential 

0.0 

3.4 

Agricultural 

20.6 

17.8 

Pasture 

5.2 

4  .2 

Institutional 

0.0 

1.8 

Roads 

0.0 

2.5 

Natural 

63.9 

60.3 

4 

10 

5946 

Low  Residential 

0.5 

0.5 

Medium  Residential 

7.2 

9.7 

Commercial 

0.7 

0.8 

Industrial 

0.4 

0.4 

Agricultural 

36.9 

34.7 

Pasture 

1.7 

1  .7 

High  Connercial 

0.0 

0.1 

Institutional 

1.3 

1.6 

Roads 

0.0 

0.4 

Natural 

51.3 

50.1 

4 

11 

3460 

Low  Residential 

1.5 

1.5 

Medium  Residential 

0.1 

0.6 

Commercial 

0.4 

0.4 

Industrial 

0.8 

2.3 

Agricultural 

10.7 

9.8 

Pasture 

13.2 

12.9 

Institutional 

0.1 

1  .0 

Natural 

73.2 

71.4 

4 

13 

3428 

Low-Medium  Residential 

20.6 

32.1 

High  Residential 

1.3 

2.7 

Commercial 

0.2 

0.5 

Industrial 

0.0 

2.8 

Agricultural 

12.4 

8.5 

Pasture 

8.3 

7.0 

High  Conmerclal 

0.0 

0.2 

Institutional 

1.5 

2.3 

Natural 

55.7 

43.1 

Open 

0.0 

1 .0 

29 


TABLE  III-l 


REACH 

SUB-BASIN 

NO. 

NO. 

LAND  USE 


Low  Residential 

Medium  Residential 

Conmerclal 

Industrial 

Agricultural 

Pasture 

Natural 


Low  Residential 

Medium  Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Institutional 

Natural 


EXIST 


11259 


Low-Medium  Residential 

High  Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Institutional 

Natural 

Open 


Low  Residential 

High  Residential 

Commercial 

Industrial 

Pasture 

Natural 


Low  Residential 

High  Residential 

Commercial 

Industrial 

Agricultural 

Pasture 

Natural 


0.0 

3.6 

30.4 

27.9 

5.7 

5.1 

0.1 

0.9 

56.6 

51.8 

0.4 

0.4 

0.4 

0.4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

48.7 

50.9 

50.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.5 

0.5 

83.1 

83.1 

16.4 

16.4 

TABLE  II 1-1 


LAND  USE  BY  WATERSHED  AND  REACH 


REACH 

SUB-BASIN 

AREA 

LAND  USE 

EXIST 

1990C 

NO. 

NO. 

(ACRES) 

% 

5 

22 

39550 

Low  Residential 

1.0 

1.0 

l 

|  Medium  Residential 

0.2 

0.2 

Commercial 

0.0 

0.0 

Industrial 

0.0 

0.0 

1 

Agricultural 

0.1 

0.1 

Pasture 

23.9 

23.9 

Natural 

74.8 

74.8 

5 

25 

9880 

Low  Residential 

0.0 

0.0 

High  Residential 

0.0 

0.0 

' 

Commercial 

0.0 

0.0 

Industrial 

0.0 

0.0 

Pasture  1 

15.1 

15.1 

Natural 

| 

84.9 

84.9 

5 

24 

1320 

Low  Residential 

0.0 

0.0 

High  Residential 

0.0 

0.0 

| 

Commercial 

0.0 

0.0 

1 

Industrial 

0.0 

0.0 

Pasture 

86.7 

86.7 

Natural 

13.3 

13.3 

RIVER  GEOMETRY 


Cross  section  data  at  irreqular  Intervals  along  the  entire  stream  system 
were  provided  by  the  Savannah  District.  The  data  were  provided  in  a  format 
for  input  to  computer  program  HEC-2,  Water  Surface  Profiles  [3].  HEC-2 
output  provided  Information  on  energy  grade  line  elevations  which  is  a 
required  Input  to  WQRRS.  More  importantly,  these  same  cross  sections  are 
Input  to  computer  program  GEDA,  Geometric  Elements  from  Cross  Section 
Coordinates  [4],  GEDA  provides  output  of  vertically  layered  geometric 
data  (l.e.,  cross  section  area,  top  width,  hydraulic  radius,  composite 
Manning's  n,  etc.)  at  regularly  spaced  nodal  points  (e.g.,  one  half  mile 
apart),  as  required  by  WQRRS. 

The  preparation  of  geometric  data  for  the  WQRRS  model  is  relatively 
automatic  once  the  basic  data  of  station-elevation  coordinate  points 
have  been  obtained  either  from  field  surveys  or  from  contour  maps. 


HYDROLOGY 


The  Oconee  basin  was  found  to  have  little  available  hydrologic  data. 

For  the  desired  watershed  modeling  purposes,  hydrologic  data  can  be 
considered  virtually  non-existent.  This  serious  lack  of  data  required 
numerous  assumptions.  The  accuracy  of  these  assumptions  cannot  be 
evaluated  except  that  the  results  did  not  seem  unreasonable  in  terms 
of  general  hydrologic  engineering  judgment. 

Only  one  USGS  stage  gage  with  hourly  flow  records  was  in  operation  during 
the  selected  study  period,  1570.  This  gage  is  located  at  river  mile  22.3 
on  the  Middle  Oconee  River.  To  obtain  the  inflow  across  the  study 
boundaries  on  the  Middle  and  North  Oconee  Rivers  (i.e.,  river  mile  31.5 
and  34.8  respectively),  the  hourly  flow  rate  at  the  USGS  gage  was  multi¬ 
plied  by  the  ratio  of  drainage  area  above  the  gage  to  that  above  each 
boundary. 

Modified  Puls  routing  criteria  for  the  North,  Middle  and  Main  Oconee 
Rivers  were  provided  by  the  Savannah  District  for  selected  control 
points.  These  data  were  linearly  interpolated  to  obtain  criteria  at 
each  load  point  (i.e.,  tributary  inflows,  sewage  treatment  plant 
effluents,  and  withdrawal  locations).  These  data  are  shown  in  Table 
1 1 1 -2 . 
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TABLE  II 1-2 

MODIFIED  PULS  ROUTING  CRITERIA 


NORTH  OCONEE  RIVER 


Bate-' 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

RIVER 

MILE 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

33.5-34.8 

16 

180 

23.2-23.7 

136 

5282 

25 

360 

253 

10600 

112 

1800 

| 

289 

13921 

31.8-33.5 

27 

180 

21 .7-23.2 

206 

5282 

42 

360 

340 

10600 

247 

1800 

425 

13921 

31.0-31.8 

14 

l 

180 

20.3-21.7 

65 

5282 

22 

360 

125 

10600 

119 

1800 

163 

13921 

28.7-31.0 

35 

180 

19.9-20.3 

35 

5312 

56 

360 

67 

10900 

352 

1800 

85 

14338 

26.0-28.7 

55 

180 

19.0-19.9 

63 

5312 

81 

360 

123 

10900 

246 

1800 

155 

14338 

24.3-26.0 

27 

180 

15.6-19.0 

440 

5312 

42 

360 

786 

10900 

174 

1800 

995 

14338 

23.8-24.3 

136 

180 

12.3-15.6 

435 

5312 

253 

360 

804 

10900 

289 

1800 

1039 

14338 

23.7-23.8 

30 

5282 

57 

10600 

65 

13921 
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TABLE  1 1 1- 2 

MODIFIED  PULS  ROUTING  CRITERIA 


MIDDLE  AND  MAIN  OCONEE  RIVER 


BH 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

31.2-31 .5 

5 

270 

16.5-17.0 

52 

7100 

7 

525 

100 

13500 

17 

2700 

122 

17000 

29.5-31.2 

521 

7100 

15.1-16.5 

407 

7100 

829 

13500 

785 

13500 

964 

17000 

961 

17000 

25.9-29.5 

1148 

7100 

12.9-15.1 

501 

7100 

2397 

13500 

961 

13500 

2937 

17000 

1391 

17000 

24.6-25.9 

145 

7100 

12.3-12.0 

50 

7100 

279 

13500 

95 

13500 

361 

17000 

120 

17000 

22.0-24.6 

380 

7100 

11.4-12.3 

353 

8675 

636 

13500 

654 

19300 

784 

17000 

819 

2  5222 

19.9-22.0 

395 

7100 

9.8-11.4 

564 

8675 

731 

13500 

i 

! 

1173 

19300 

885 

17000 

i 

1462 

25222 

17.3-19.9 

571 

7100 

8. 9-9. 8 

214 

8675 

1034 

13500 

364 

19300 

1247 

17000 

433 

25222 

17.0-17.3 

52 

7100 

8. 0-8. 9 

133 

8675 

100 

13500 

206 

19300 

122 

17000 

235 

25222 

TABLE  1 1 1-2 

MODIFIED  PULS  ROUTING  CRITERIA 


MIDDLE  AND  MAIN  OCONEE  RIVER 


STORAGE 

(AF) 

OUTFLOW 

(cfs) 

RIVER 

MILE 

STORAGE 

(AF) 

OUTFLOW 

(cfs) 

.0 

398 

8675 

0.8-2. 8 

359 

8675 

777 

19300 

1031 

19300 

957 

25222 

1249 

25222 

.2 

734 

8675 

0.0-0. 8 

360 

8675 

1318 

19300 

983 

19300 

1630 

25222 

1248 

25222 

.8 

571 

8675 

1440 

19300 

1742 

25222 

WATER  QUALITY 


The  Oconee  basin  was  found  to  have  little  available  water  quality  data  of 
practical  use  to  this  study.  It  was  originally  thought  that  adequate  water 
quality  data  were  available  to  calibrate  the  water  quality  simulation  models, 
STORM  and  WQRRS.  This  serious  lack  of  data  required  numerous  assumptions  which 
could  not  be  verified  with  field  data.  The  only  evaluation  that  could  be  made 
was  that  the  results  appeared  reasonable  In  the  light  of  other  experience. 

Since  1968,  water  quality  data  have  been  collected  at  the  USGS  stage  gage 
on  the  Middle  Oconee  River  (l.e.,  river  mile  22.3)  and  since  1974  at  the 
Athens  water  Intake  on  the  North  Oconee  River.  The  data  from  these  two 
sampling  locations  together  with  the  data  from  Smith  [5],  Appendix  A,  were 
used  to  estimate  base  flow  quality  data  for  the  boundary  condition  and  the 
tributaries.  The  boundary  quality  condition  was  held  constant  during 
storm  events  and  the  tributary  Inflow  quantity  and  quality  was  obtained 
from  the  STORM  model  output. 


IV  MODELING  CONCEPTS  APPLIED 


STORM 

The  Storage,  Treatment,  Overflow  Runoff  Model  (STORM)  Is  a  continuous 
simulation  model  designed  to  be  used  in  metropolitan  master  planning 
studies  for  evaluating  storage  and  treatment  capacities  required  to 
reduce  overflows.  Pollutograph  (pollutant  mass-emmisslon  rates)  loadings 
can  also  be  computed  for  use  in  a  receiving  water  assessment  model. 

Since  STORM  is  intended  for  use  in  planning  studies  or  for  screening 
alternatives,  some  of  its  analytical  techniques  are  necessarily 
simplified.  For  example,  the  two  procedures  used  to  compute  the  quantity 
of  runoff  are  the  coefficient  method  and  the  United  States  Soil  Conservation 
Service  (SCS)  method.  In  the  coefficient  method,  a  single  land-use  weighted 
runoff  coefficient  is  applied  to  each  hour  of  rainfall  excess  above  depression 
storage  to  compute  runoff.  The  runoff  coefficient  is  a  function  of  only  the 
respective  runoff  coefficients  for  the  pervious  and  impervious  areas  of  the 
watershed.  Antecedent  conditions  and  rainfall  intensity  are  not  taken  Into 
account  using  this  method. 

The  SCS  runoff  curve  number  technique  is  considered  to  be  conceptually  more 
correct  than  the  coefficient  method.  The  SCS  curve  consists  of  a  nonlinear 
relationship  between  accumulated  rainfall  and  accumulated  runoff.  Since 
STORM  requires  a  continuous  analysis,  a  nrocedure  has  been  added  that 
computes  the  curve  number  for  each  event  based  on  the  number  of  dry  hours 
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since  the  previous  runoff  event  and  prior  evapotranspiratlon  and  percolation. 
Unit  hydrographs  can  be  used  to  transform  the  surface  runoff  excesses  Into 
basin  outflow  hydrographs. 


Loads  and  concentrations  for  six  basic  water  quality  parameters  are  computed. 
These  are  suspended  and  settleable  solids,  biochemical  oxygen  demand,  total 
nitrogen,  total  orthophosphate,  and  total  col i forms.  Urban  and  nonurban 
areas  may  be  described  by  up  to  20  land  uses.  Other  features  of  STORM  are 
the  capabilities  to  compute  snowfall/snowmelt,  dry-weather  flow  quantity 
and  quality,  and  land  surface  erosion. 

STORM  has  a  unique  advantage  of  being  able  to  accept  discharge  hydrographs 
as  input  for  computing  the  associated  wash-off  of  constituents. 

WQRRS 

The  Water  Quality  for  River-Reservoir  Systems  (WQRRS)  model  has 
capability  for  ecologic  evaluation  of  rivers  or  reservoirs.  It  is  a 
dynamic  continuous  simulation  model.  The  model  consists  of  three  separate 
but  integrable  modules.  These  are  the  reservoir  module,  the  stream 
hydraulic  module,  and  the  stream  quality  module.  Since  each  module  is  a 
stand-alone  proqram,  the  reservoir,  the  streamflow  routing,  or  the  stream 
water  quality  module  may  be  executed,  analyzed  and  interpreted  independently. 
The  three  computer  proorams  may  also  be  integrated  into  a  complete  river 
basin  water  quality  analysis. 

The  reservoir  section  of  the  proqram  estimates  the  water  quality  condition 
in  deep  impoundments  that  can  he  represented  as  one-dimensional  systems 
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In  which  the  Isotherms,  or  contours  of  any  parameter,  are  horizontal. 

This  approximation  Is  generally  satisfactory  in  lakes  with  lonq  residence 
times.  However,  the  approximation  is  less  satisfactory  in  shallow 
Impoundments  or  those  that  have  a  rapid  flow-through  time.  Systems 
that  have  a  rapid  flow-through  time  are  often  fully  mixed  and  can  be 
treated  as  slowly  moving  streams  using  the  stream  section  of  the  model. 

The  stream  hydraulic  section  of  the  model  includes  six  hydraulic  calcu¬ 
lation  options.  This  module  Is  capable  of  handling  hydraulic  behavior 
for  both  the  "gradually  varied"  steady  and  unsteady  flow  regimes.  Peak 
flows  from  storm  water  runoff  or  irregular  hydropower  releases  can  be 
accuratel y  represented . 

In  the  stream  quality  module,  the  rate  of  transport  of  quality  parameters 
can  be  accurately  represented  and  peak  pollutant  loads  into  the  steady  or 
unsteady  hydraulic  environment  can  be  simulated.  The  stream  portions  of 
WQRRS  have  two  automatic  interface  options  for  use  with  the  STORM  model. 


V  MODELING  RESULTS 


STORM  RUNOFF  QUANTITY  AND  QUALITY 

The  approach  used  to  calculate  storm  runoff  was  to  sub-divide  the  total 
study  area  Into  a  number  of  Individual  watersheds  and  apply  STORM  to  each. 
Criteria  affecting  the  number  of  watersheds  Include  the  degree  of  refinement 
In  discrete  points  along  the  receiving  water  body  where  Individual 
calculations  are  to  be  made  and  manageability  of  data  for  the  entire 
study.  A  total  of  32  Individual  watersheds  were  Identified  and  are  shown 
on  Figure  III-2.  Several  of  the  watersheds  In  the  Savannah  District  Data 
Bank  were  further  subdivided  so  as  to  provide  better  definition  of  quality 
In  the  urbanized  river  reaches. 

The  first  major  effort  in  the  STORM  application  was  to  assemble  and  edit 
hourly  precipitation  data.  Five  recording  rain  gages  exist  In  or  near 
the  study  area;  the  locations  are  shown  In  Figure  lll-Z.  Only  the  Athens 
gage  (Station  No.  0435)  was  a  Class  I  U.S.  Weather  Service  gage.  The 
other  four  gages  are  supplemental  locations  and,  as  a  result,  the  data 
have  not  been  corrected  for  gage  failures.  Numerous  places  existed  on 
the  tapes  where  the  gage  had  failed  and  the  accumulated  precipitation  was 
shown  In  the  first  hour  of  resumption  of  recording.  A  special  editor 
program  was  written  to  locate  these  gage  failures  and  redistribute  the 
precipitation  to  the  hours  in  which  they  occurred.  The  precipitation 
was  distributed  evenly  over  the  hours  of  gage  failure. 


The  continuous  precipitation  histories  (1948-72)  were  used  to  assess  the  average 
annual  land  surface  erosion  for  both  existing  and  future  conditions.  The  single 
year  1970  was  used  for  storm  water  simulations  since  it  was  the  month  of 
October  1970  that  was  studied  in  the  receiving  water  analysis  using  WQRRS. 

No  data  existed  with  which  to  calibrate  the  rainfall -runoff  calculations  in 


STORM  for  the  Oconee  study.  The  various  soil  moisture  characteristics  required 

for  the  SCS  runoff  method  were  estimated.  The  October  percentage  runoff  for 

several  nearby  gages  served  as  a  guide.  The  tributary  flows,  when  combined 

and  routed  to  a  gage  near  the  downstream  study  boundary  showed  fair  agreement 

with  the  observed,  however  no  check  could  be  made  on  Individual  tributary 

flows.  The  estimated  runoff  characteristics  are  shown  below: 

Soil  Moisture  Max  Initial 

at  Saturation  Abstraction  Capacity 

LAND  USE  (SMAX),  Inches  (DEPR),  Inches 

Open 

14.40 

2.88 

Low  D.  Residential 

11.70 

1.17 

Medium  D.  Residential 

7.54 

0.11 

Low-Medium  D.  Residential 

9.62 

0.18 

High  D.  Residential 

4.28 

0.06 

Agricultural 

7.54 

1.51 

Industrial 

3.33 

0,05 

Commercial 

2.34 

0.04 

Pasture 

10.00 

2.00 

High  D.  Commercial 

1.24 

0.02 

Institutional 

3.16 

0.05 

Roads 

4.70 

0.07 

Natural 
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No  data  existed  with  which  to  calibrate  the  runoff  quality  calculations  In 
STORM  for  the  Oconee  Study.  The  various  pollutant  accumulation  rates  required 
to  regulate  the  quality  were  estimated  based  on  data  from  the  literature 
[6,7].  Adjustments  were  made  during  calibration  so  that  tributary  storm 
water  concentrations  for  existing  conditions  did  not  greatly  exceed  certain 
measured  concentrations  in  the  river.  Storm  water  quality  calculations  were 
not  calibrated  directly  to  the  river  concentration  for  three  reasons: 

1)  minimal  data  existed,  2)  the  data  consisted  of  grab  samples  taken  at 
Infrequent  Intervals,  and  3)  there  were  no  indications  that  the  measurements 
were  taken  during  periods  of  tributary  storm  runoff.  Table  V-l  shows  the 
adopted  pollutant  accumulation  rates. 


TABLE  V-l 

Pollutant  Accumulation  Rates  (lb/ac/day) 


Land 

Use 

Susp 

Solids 

11 

mum 

B0D5 

N 

— 

P04 

Col i form 
109MPN/ac/day 

Low  Res 

.12 

.09 

.002 

.0002 

.0004 

.60 

LM  Res 

.43 

.16 

.004 

.0008 

.0006 

.62 

Med  Res 

.45 

.18 

.004 

.0008 

.0006 

.63 

High  Res 

3.10 

.99 

.006 

.0006 

.0020 

4.9 

Coml 

3.60 

1 .35 

.022 

.0060 

.0040 

4.5 

Ind 

6.00 

2.25 

.020 

.0055 

.0030 

5.0 

Agr 

7.20 

2.70 

.001 

.0012 

.00002 

.25 

Pasture 

.24 

.10 

.001 

.0002 

.0002 

.50 

HI  Coml 

3.90 

1.44 

.016 

.0065 

.0048 

5.0 

Instl 

3.10 

1.17 

.006 

.0006 

.0020 

6.0 

Natural 

.10 

.04 

.001 

.0001 

.000002 

.0005 

Open 

.24 

.08 

.001 

.0002 

.0002 

.50 
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Dry  weather  sewage  flow  was  simulated  for  those  basins  with  significant 
urban  land  use.  Dry  weather  flow  option  three  was  used  since  it  allows 
computations  to  be  made  on  the  basis  of  land  use  and  population.  The 
coefficients  used  are  shown  in  Table  V-2.  Domestic,  Commercial  and 
Industrial  coefficients  were  taken  from  References  8  and  9,  with  some 
minor  modifications.  Pipe  Infiltration  coefficients  were  estimated  so 
that  the  quality  concentrations  did  not  exceed  those  for  baseflow  from 
non-urban  subbasins.  The  coefficients  were  assumed  to  remain  constant 
for  the  alternative  future. 


TABLE  V-2 

Dry-Weather  Flow  Coefficients 


Domestic 

Commercial 

Industrial 

Infiltration 

low 

|mgd/acre) 

lOO^ 

.03 

.01 

.0005 

Suspended  Solids 
(lb/day/ac) 

1.3^ 

1.9 

2.6 

.08 

iettleable  Solids 
(lb/day/ac) 

.22-/ 

.33 

.44 

.008 

10°5 

(lb/day/ac) 

.20^ 

.30 

.40 

.002 

1 

(lb/day/ac) 

.04—^ 

.05 

.06 

.0012 

*4 

(lb/day/ac) 

.oi  y 

.012 

.02 

.0004 

Collform 

(lO^MPN/day/ac) 

.64^/ 

.86 

.86 

.0125 

1/  gallons/day/caplta  for  Domestic 
2J  pounds/day/capita  for  Domestic 
3/  109  MPN/day/capIta  for  Domestic 


Since  the  dry  weather  flow  algorithm  In  STORM  calculates  loads  and  con¬ 
centrations  of  raw  waste  water,  reductions  must  be  made  to  account  for 
treatment  that  exists  In  the  study  area.  An  assumption  was  made  that  the 
level  of  treatment  remains  constant  for  the  alternative  future.  The 
following  removal  efficiencies  were  used  for  each  subbasin  having  dry  weather 
flow. 


Treatment  Efficiencies  Used  in  STORM 
(percent) 


Suspended  Solids 

87 

bod5 

87 

Orthophosphate 

80 

Settleable  Solids 

87 

Nitrogen 

80 

Collform 

92 

Table  V-3  shows  predicted  tributary  loads  in  pounds  for  the  10  month 
period  of  January  1970  through  October  1970.  While  these  loads  cannot 
be  used  as  an  evaluation  objective  In  themselves,  they  are  useful  to 
distinguish  trends.  In  every  case  the  predicted  loads  for  1990C  land  use 
pattern  exceeded  those  for  existing  conditions.  These  loads  were  not  used 
for  the  Instream  analysis.  The  receiving  water  analysis  was  accomplished 
using  hourly  loads  and  concentrations  for  the  month  of  October  1970. 

Storm  water  guantity  and  guality  were  also  simulated  for  the  Pendergrass 
detailed  study  area,  however  since  a  receiving  water  analysis  was  not 
performed  there  was  no  need  to  predict  individual  subbasin  loadings. 

Table  V-4  summarizes  the  predicted  storm  runoff  guality  loadings  for  the 
Pendergrass  detailed  study  area. 
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TABLE  V-3 

Predicted  Washoff  Of  Pollutants 
January  through  October  1970 
For  Athens,  Georgia 


Sub-  Land  Use 
Basin  Condition 


Suspended  Settleable  BODc  N  PO*  Coll  form 

(lbs)  (lbs)  (lbs)  (lbs)  (lbs)  (109  MPN) 


TABLE  V-3  (Cont) 


Predicted  Washoff  Of  Pollutants 
January  through  October  1970 
For  Athens,  6eorq1a 


Sub 

Basin 

Land  Use 
Condition 

Suspended  Settleable 
(lbs)  (lbs) 

BOD, 

(lbf) 

N 

(lbs) 

P04 

(lbs) 

Collform 
(109  MPN) 

10 

Exist 

289,780 

37,070 

17,960 

4,613 

1,344 

172,652 

1990C 

333,103 

42,923 

20,516 

5,268 

1,565 

246,885 

11 

Exist 

32,252 

3,586 

1,516 

364 

102 

30,428 

1990C 

41 ,906 

5,182 

5,482 

604 

185 

84,808 

12A 

Exist 

117,069 

14,402 

7,755 

1,898 

605 

226,020 

1990C 

207,271 

25,616 

13,865 

3,401 

1,161 

504,939 

12B 

Exist 

751 ,269 

145,068 

49,750 

12,276 

3,920 

1,089,028 

1990C 

735,841 

144,459 

48,679 

11,989 

3,857 

1,144,164 

12C 

Exist 

229,663 

30,845 

17,105 

4,239 

1,216 

439,079 

1990C 

686,401 

102,465 

48,927 

12,159 

4,272 

1,001,348 

13 

Exist 

106,697 

14,890 

9,531 

2,934 

445 

175,192 

1990C 

200,750 

27,880 

19,367 

5,209 

1,792 

458,203 

14 

Exist 

237,387 

29,749 

14,208 

3,649 

1,060 

128,786 

1990C 

251 ,051 

31,778 

15,142 

3,893 

1,143 

149,505 

15 

Exist 

16,468 

1 ,982 

873 

197 

52 

14,138 

1990C 

23,865 

2,950 

1,291 

297 

85 

27,180 

16 

Exist 

235,374 

30 ,904 

13,910 

4,540 

928 

171,134 

1990C 

508,469 

64,824 

33,360 

9,468 

2,571 

872,182 

TABLE  V-4 


Storm  Runoff  Quality  Loadings 
Pendergrass  Study  Area  Jan  -  Oct  1970 


Suspended 
Solid s  (lb) 

Settleable 
Sol ids(lb] 

b°d5 

Ob? 

Nitrogen 

Ob) 

PO4 

Ob) 

Col  1  form 
(109MPN) 

Existing 

195,244 

23,730 

12,808 

3,099 

933 

214,766 

1990  B 

1,113,700 

142,800 

74,755 

18,467 

6,234 

2,139,334 
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LAND  SURFACE  EROSION 


Land  surface  erosion  yield  was  computed  by  the  Universal  Soil  Loss  Equation. 

The  equation,  as  implemented  in  STORM,  requires  a  continuous  hourly  precipitation 
record  to  serve  as  the  prime  mover  In  the  analysis.  The  period  of  January 
1948-December  1972  at  the  Athens  gage  for  the  Athens  area  (Winder  gage  for 
Pendergrass)  was  used  for  both  existing  and  future  conditions.  The  K,  LS,  C, 

P  and  SDR  terms  in  the  equation  are  shown  in  Table  V-5  [10,  11]. 

The  average  annual  land  surface  erosion  was  computed  for  the  period  of  record 
and  from  several  trial  runs  It  was  determined  that,  for  the  period  of  January 
1,  1961  to  December  31,  1961  the  land  surface  erosion  approximated  the  average 
annual  land  surface  erosion  for  the  period  of  record.  This  shorter  period 
was  used  in  subsequent  runs  to  calculate  the  average  annual  land  surface 
erosion  for  the  various  subbasins  in  the  Athens  Study  Area. 

An  Important  consideration  in  the  land  surface  erosion  analysis  was  the 
effect  of  exposed  soil  in  areas  under  development.  For  each  grid  cell  the 
land  use  for  existing  and  future  conditions  were  compared  and  the  number  of 
cells  with  changed  land  use  were  counted.  It  was  then  assumed  that  the 
change  In  land  use  will  be  uniformly  distributed  over  the  15  year  period 
1976-1990.  Therefore,  the  area  under  development  for  any  one  year  is 
approximately  7%  of  the  total  change  during  the  15  years.  For  that  area 
under  development  the  factors  representing  the  soil  cover  were  modified 
to  reflect  denuded  soil.  Specifically,  the  Cover  Factor  and  Erosion 
Control  Factor  were  set  to  100. 
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TABLE  V-5 


Soil  Erodibility  Factors 


SOIL 

NO. 

SOIL 

CODE 

SERIES 

NAME 

K 

SOIL 

NO. 

■a 

SERIES 

NAME 

a 

2 

Ak 

39 

Ln 

Loulsburg 

0.24 

3 

Am 

41 

Me 

4 

An 

Appling 

0.32 

43 

Mg 

Madison 

0.32 

5 

As 

44 

Mi 

Madison 

0.32 

6 

Ax 

Appling 

0.32 

45 

Mm 

Madison 

0.32 

7 

Bfs 

Buncombe 

0.17 

Mm 

Louisa 

0.28 

8 

Ca 

47 

Mv 

Musel la 

0.28 

9 

Cb 

Cecil 

0.32 

48 

Pa 

10 

Ce 

49 

Pf 

Pacolet 

0.32 

11 

Cf 

50 

Pg 

Pacolet 

0.32 

12 

Ci 

Colfax 

51 

Ph 

Pacolet 

0.32 

13 

Coa 

Conga ree 

52 

Pi 

Pacolet 

0.32 

14 

Cob 

Chewacla 

54 

Pt 

17 

Cy 

55 

Rc 

18 

CY 

Cecil 

0.32 

56 

Rok 

Rock 

0.00 

19 

CZ 

Cecil 

0.32 

57 

Tf 

21 

Oh 

Davidson 

0.32 

58 

To 

22 

Dq 

Davidson 

0.32 

61 

Wq 

23 

EW 

62 

Wk 

Worsham 

25 

Ge 

63 

Wos 

Wehadkee 

29 

Gr 

65 

LD 

Louisburq 

0.25 

31 

He 

33 

H1  1 

NOTE:  SOILS  WITH  NO  "K”  VALUE  USED  THE  DEFAULT  OF  0.32 
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TABLE  V-5 


< 


Length  for  LS  Factor 


SLOPE  * 


LENGTH  (ft) 


0  -  2.00  200 

2.01  -  6.00  275 

6.01  -  10.00  175 

10.01  -  15.00  75 

15.01  -  25.00  50 


Cover  and  Erosion  Control  Factors 


Land 

Use 

Cover  Factor  (C) 

Erosion  Control 
Factor  %  (P) 

1. 

Open 

1.3 

05 

2. 

Low  Residential 

0.3 

85 

3. 

Medium  Residential 

0.3 

70 

4. 

High  Residential 

0.3 

60 

c 

J  • 

Agricultural 

40.0 

95 

6. 

Industrial 

10.0 

40 

7. 

Commercial 

1.2 

20 

8. 

Pasture 

2.0 

95 

9. 

High  Commercial 

1.0 

10 

10. 

Institutional 

10.0 

40 

11. 

Roads 

5.0 

60 

12. 

Natural 

0.3 

95 

TABLE  V-5 


Sediment  Delivery  Ratio  (SDR) 


WATERSHED  NO.* 

AREA,  ac 

SDR 

1 

6669 

0.18 

2 

10444 

0.16 

3 

2272 

0.23 

4 

7915 

0.17 

5 

7648 

0.17 

6 

16301 

0.14 

7 

6101 

0.18 

8 

12164 

0.16 

9 

6242 

0.18 

10 

1290 

0.25 

11 

5946 

0.19 

12 

3460 

0.21 

13 

3438 

0.21 

14 

5039 

0.19 

15 

3103 

0.21 

16 

11253 

0.16 

Pendergrass 

7067 

0.17 

All  other  factors  In  the  soil  loss  equation  remained  the  same  as  In 
the  developed  condition.  The  predicted  land  surface  erosion  for  the 
Athens  and  Pendergrass  study  areas  are  shown  In  Tables  V-6  and  V-7. 
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Table  V- 6 


Athens*  Georgia 


Average  Annual 

Land  Surface  Erosion  (tons) 

Watershed  No 

.*  Existing 

1990  C 

1 

76900 

72000 

2 

112000 

115300 

3 

5000 

6300 

4 

85000 

83900 

5 

51100 

55700 

6 

118700 

306200 

7 

10100 

16900 

8 

222700 

219900 

9 

93100 

88500 

10 

17900 

17100 

11 

96800 

93800 

12 

24500 

24400 

13 

23500 

23100 

14 

89000 

67000 

15 

4800 

5200 

16 

163100 

157400 

♦Savannah  District  Watershed  Identification 
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TABLE  V-7 


Average  Annual 
Land  Surface  Erosion 
Pendergrass  Area 


Existing  Land  Use  105220  tons 

Alternative  B  Land  Use  94910  tons 
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RECEIVING  WATER 


Analysis  of  Existing  Condition 

The  WQRRS  model  accepts  Input  tributary  conditions  derived  using  the  STORM 
model  on  each  of  the  32  subbasins  (see  Figures  III-2,  III-3  and  III-4) 
draining  Into  the  portion  of  the  Oconee  River  within  the  specified  study 
boundaries  and  Imposes  these  loadings  on  a  base  flow  condition.  The  two 
sewage  treatment  plants  and  the  Athens  water  Intake  are  accounted  for  based 
on  mean  monthly  data  from  the  State  of  Georgia,  except  for  unmeasured 
parameters  which  were  then  estimated  from  textbook  average  conditions  [12]. 
These  Input  data  are  shown  In  Table  V-8.  An  accounting  Is  made  of  the  mass 
balance  at  each  tributary  location  and  the  resulting  mixture  Is  transferred 
(l.e.,  routed)  downstream  to  the  next  tributary  location  with  the  proper 
reactions  and  Interactions  being  calculated  according  to  the  estimated 
travel  time  between  tributaries  and  the  Input  system  coefficients  shown 
In  Table  V-9. 

The  Initial  quality  condition  for  selected  locations  Is  shown  In  Table  V-10. 
The  quality  at  all  other  locations  Is  obtained  by  linear  Interpolation.  The 
values  shown  for  river  mile  34.8  and  31.5  on  the  North  Oconee  and  Middle 
Oconee  respectively  are  also  the  base  flow  quality  conditions  which  enter 
the  study  area  at  the  upper  boundaries. 

The  base  flow  condition  on  the  tributaries  during  non-storm  periods  Is 
dependent  on  the  proportion  of  the  drainage  area  having  residental  land 
use.  Table  V-U  shows  the  tributary  base  flow  used  as  Inflow  to  the  Oconee 
during  non-storm  periods. 
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TABLE  V-8 

INPUT  DATA  FOR  SEWAGE  TREATMENT  PLANTS  AND 
WATER  TREATMENT  PLANTS  FOR  EXISTING  AND 
ALTERNATIVE  FUTURE  C  LAND  USE 


Parameter  1/ 

North  Oconee 
STP(R.M.  19.9) 

Exist  "  Alt.C 

Middle  Oconee 
STPfR.M.  17.3) 
Exist  JCTt.C 

Athens  Water 
Intake(R.M.  23.7) 
Exist  Alt.C 

Q  (cfs)  y 

10.1 

15.3 

3.25 

5.0 

16.6 

24.9 

Temperature  (°C) 

3 / 

3/ 

3/ 

If 

DO  (assumed) 

0 

0 

0 

0 

BODS  4/ 

99 

99 

84 

84 

Collform  (assumed) 
(MPN/100  ml) 

ZOO 

200 

200 

200 

Detritus 

(25%  of  susp.  solids) 

8.25 

8.25 

7.25 

7.25 

NH3  [5] 

10 

10 

10 

10 

N03  [5] 

20 

20 

20 

20 

N02  [5] 

.05 

.05 

.05 

.05 

P04  [5] 

12 

12 

12 

12 

TDS 

244 

244 

154 

154 

Algae  (assumed) 

.001 

.001 

.001 

.001 

Zooplankton 

(assumed) 

.001 

.001 

.001 

.001 

pH  (units) 

6.7 

6.7 

7.3 

7.3 

Akallnlty  (assumed) 

100 

100 

100 

100 

1/  mg/1  except  as  noted. 

y  Flow  for  Alternative  C  equals  existing  flow  times  estimated. proportional  Increase 
In  population. 

y  Water  temperature  equals  mean  dally  air  temperature  minus  2*C  [13]  except  during 
storm  events  when  water  temperature  equals  the  hourly  air  temperature. 

4/  Uncorrected  for  NHj,  N02  and  Detritus  oxygen  demand  Cl 3 3 • 
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TABLE  V-9 

INPUT  SYSTEM  COEFFICIENTS 

REACTION  PATF  MULTIPLIER  PARAMETERS 


CALIBRATION  MAGNITUDES  CALIBRATION  TEMPERATURES 


Kl 

K2 

K3 

K4 

T  1 

T  2 

T  3 

T  4 

ALGAE  1 

.  10 

.98 

.98 

.  10 

5.0 

22.0 

25.0 

34.0 

ALGAE  2 

.10 

.98 

«  9  R 

.10 

10.0 

28.0 

30.0 

40.0 

ZOOPLANKTON 

.10 

.98 

.98 

.10 

5.0 

28.0 

30.0 

39.0 

BENTHIC  ANIMALS 

.  10 

.9e 

.98 

.  10 

5.0 

22.0 

25.0 

33.0 

FISH  1 

.10 

.98 

.98 

.  10 

5.0 

20.0 

20.0 

25.0 

FISH  2 

.10 

.9  8 

.98 

.10 

10.0 

27.0 

30.0 

33.0 

FISH  3 

.  10 

.98 

.98 

.10 

5.0 

22.0 

30.0 

36.0 

BOD 

.10 

.98 

4.0 

30.0 

NH  3-N 

.10 

.98 

4.0 

30.0 

N02-N 

.10 

.98 

4.0 

30.0 

DETRITUS 

.10 

.98 

4.0 

30.0 

DECAY  CCJEFf  ICIENTS,  PER  DAY 

MAX  VALUE 

BOD 

.  100 

NH3-N 

.050 

N02-N 

.  2C0 

DETRITUS 

.001 

COL  I  FOP  M  (AT  20  OEG  C) 

.500 

010  TEMPERATURE  COEFFICIENT  FCR  CCLIFCRM  1.040 

CHEMICAL  COMPOSITIONS  OF  BIOTA 

C  N  P 


ALGAE 

.  500 

.090 

.012 

ZOOPLANKTON 

.500 

.090 

.012 

FISH 

.500 

.090 

.012 

BENTHOS 

.500 

.090 

.012 

detritus 

.500 

.090 

.01 2 

TABLE  V-9  (cont'd) 


RESPIRATION  PATES, PFP  DAY 

MAX  VALUE 

PHYTOPLANKTON 

. 500E-01 

ZOOPLANKTON 

. 200E-01 

FISH 

. I00E-02 

BENTHOS 

. 10CF-02 

DETRITUS  SETTLING,  METERS/OAY 

. 1 5000 

OTHER  PHYTOPLANKTON  DATA 

SETTLING,  METER/DAY 

.15000 

.15000 

OXYGENATION  FACTOR 

1.600 

PRFFERANCE 

.670 

.330 

SELFSHADING  PEP  “G/L/M 

0 

MAXIMUM  SPECIFIC  GROWTH  PATE, 

PER  DAY 

PHYTOPLANKTON,  2  GROUPS 

.  1  COF  0 1 

•  200F+-01 

ZOOPLANKTON 

. 15CE+00 

FISH,  3  GROUPS 

•200E-01 

.250E-01 

.  200E-0 1 

BENTHOS 

. 2C0F-01 

HALF-SATURATION  CONSTANTS  CF 

ALGAC 

LIGHT 

Cn2 

N 

PJ4 

ALGAE  L 

.003 

.  020 

.200 

.  030 

ALGAE  2 

.005 

.020 

.100 

.050 

HALF-SATURATIGN  CONSTANTS  for 

ZOO,  FISH  AND  BENT  HO 

ZOO  GRAZE  ON  A L G A F 

.550 

FISH  1  GPAZE  GN  ZOO 

.  05 0 

FISH  2  GRAZE  ON  ZOO 

.050 

FI SH  GRAZE  ON  RENTmCS 

500.000 

BFNTHOS  GRAZE  ON  SFD*T 

50 . C  00 

STOICHIOMETRIC  EQUIVALENCE 

02/NH3 
C2/N02 
02/0E  TR ITUS 
02/BIOMASS 
C02/B00 


OF  CHfviCAL  TRANSFORMATION 

3.300 
I  .200 
2.000 
2.000 
.200 
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TABLE  V-10 

INITIAL  QUALITY  CONDITION 


Location 

Parameter 

Magnitude  1/ 

Source 

North  Oconee-RM  34.8 

BODS 

.5 

Smith  [5] 

Detritus 

5 

Sediment 

(gm/m2) 

16 

O 

Benthos  (gm/m  ) 

.9 

NH- 

no3 

.03 

.24 

State  of  Georgia  gage 
at  Athens  Intake, 
average  1974 

NOp 

.01 

P04 

.03 

pH  (pH  units) 

7.1 

Alkalinity 

36 

Coliform 

430 

Athens  Intake 

(MPN/100  ml) 

10/24/74 

Temp.  (°C) 

f(alr)— ^ 

Willey  A  Huff  Tl3] 

DO 

8 

Assumed  at  80^  of 

D0sat  at  15°C 

TDS 

100 

Assumed 

Algae 

.001 

Zooplankton 

.001 

Fish  1  (Kg/mi) 

10 

Smith  [5] 

Fish  2  (Kg/ml) 

30 

Fish  3  (Kg/ml) 

40 

60 


TABLE  V-10 (cont'd) 
INITIAL  QUALITY  CONDITION 


>^rv. 


OP  fc**- 


TABLE  V-10  (cont'd) 
INITIAL  QUALITY  CONDITION 


Location 

Parameter 

Magnitude  1/ 

Source 

Middle  Oconee-RM  31.5 

B0D5 

.5 

Smith  [5] 

Detritus 

5 

Sediment 

(gm/m2) 

16 

p 

Benthos  (gm/m  ) 

.9 

NH, 

.10 

USGS  gage  on  Middle 

.60 

Oconee  NH3:  2/09/70 

1U3 

NO3  +  NO2:  avg  of 

NO, 

.01 

11/18/7 0  and  5/26/70 

PO 

.06 

PO4:  9/02/70  pH  and 

KU4 

Coliform:  11/18/70 

pH  (pH  units) 

7.5 

Col i form 
(MPN/100  ml) 

930 

Alkalinity 

25 

Alkalinity:  avg  of 
9/02/70  and  11/18/70 

Temp.  (°C) 

f(air)^ 

DO 

8 

80%  of  D0sat  at  1 5°C 

TDS 

100 

Assumed 

Algae 

.001 

Zooplankton 

.001 

Fish  1  (kg/ml) 

10 

Smith  [5] 

Fish  2  (kg/mi) 

30 

Fish  3  (kg/mi) 

40 

6  2 


TABLE  V-10  (cont’d) 
INITIAL  QUALITY  CONDITION 


Location 

Parameter 

Magnitude  1 J 

Source 

Middle  Oconee-RM  17.5 

B0D5 

NH- 

»°3 

*>4 

Collform 

All  other  parame 

.5 

.08 

.22 

.10 

660 

:ers  same  as  R.M.  31 . 

STORM  Base  Flow 

5 

Middle  Oconee-RM  17.0 

Detritus 

(gm/m2) 

Sediment 

(gm/m2) 

All  other  parame 

5.5 

18.4 

;ers  same  as  R.M.  17. 

Smith  [5] 

5 

Main  Oconee-RM  0.0 

All  parameters  a 

*e  same  as  Middle  Occ 

nee  R.M.  17.0 

1/  mg/1  except  as  noted 

2/  water  temperature  equals  hourly  air  temperature  during  a  storm  and  equals  mean 
daily  air  temperature  minus  2°C  during  non-storm  periods  [13]. 
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The  Initial  and  base  flow  quality  conditions  are  arbitrarily  accepted  base 
conditions  since  essentially  no  gaged  data  exists  for  the  study  period. 

All  the  final  results  must  be  Interpreted  relative  to  this  base  condition 
since  most  of  the  water  quality  calculations  are  non-linear  (l.e.,  effect  of 
saturation  values  and  temperature  corrections  on  all  reaction  rates). 

Table  V-12  shows  an  example  of  a  statistical  summary  of  the  water  quality 
condition  at  a  random  point  along  the  river  for  gristing  land  use.  The 
critical  values  (l.e.,  maximum  or  minimum)  for  some  of  these  parameters  at 
various  locations  along  the  river  have  been  plotted  in  Figures  V-l  to  V-6 
to  show  a  river  water  quality  profile  for  the  most  critical  condition 
occurring  during  the  period  simulated  for  dissolved  oxygen  (DO),  ammonia 
(NHj),  nitrate  (N03),  phosphate  (PO^),  log  coliform  bacteria  and  5-day 
carbonaceous  biochemical  oxyqen  demand  (B0D5). 

The  existing  water  quality  condition  seems  to  meet  all  the  Georgia  State 
Water  Quality  Standards  (i.e..  Table  V-l 3)  except  for  coliform  bacteria 
which  may  exceed  the  standards  5-10%  of  the  time,  from  river  mile  25  to 
12.3  respectively  on  the  North  Oconee  and  throughout  the  study  length  of 
the  Middle  and  Main  Oconee  about  10-16%  of  the  time. 

The  impact  on  the  water  quality  in  the  Oconee  due  to  the  various  tributaries 
and  sewage  treatment  plants  is  shown  in  Figures  V-l  to  V-6.  Major  point 
source  Impacts  on  the  North  Oconee  River  are  summarized  in  Table  V-14,  and 
on  the  Middle  and  Main  Oconee  River  in  Table  V-l 5. 
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Remarks  in  Tables  V-14  and  V-15  concerning  nutrients  having  significant 
Impact  refer  to  the  potential  impact  on  algae  production  in  non-turbid 
water.  Unless  significant  improvement  occurs  in  the  turbidity  of  the 
Oconee  River,  this  potential  will  be  not  be  realized. 

Sample  graphical  results  of  the  simulations  are  shown  in  Appendix  B. 
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TABLE  V-12 

WATER  QUALITY  AT  RIVER  MILE  11.5 
EXISTING  UNO  USE 


POST-PHOCESSuW  FUR  -aw RS  APRIL  1976 

MYOROLOG IC  EmGINFFHIMG  CENTER  DAVIS,  CA 


ucgnee  river  water  'Jijalitv  study  *»wqr«s  statistical  post-processor** 

REACH  a  MIOOLE  UCUNEf  RIVER  (R.m.  17,0-8.01 _ 

QUALITY  data  RaSLC'  Ur.  ExISTTnG  U Aiwu  USE 

•  ***•«*•*•*•**«*•••  I»,PUT  DATA  MtMtMUAMMtttMM 


BEGINNING  OF  REACH  RtvEL  n  tLE  17,00 

END  OF  REACH  RIVER  MILE  _ 8,00  _ 

SURREACH  LENGTH  (MILES)  ,50 

COMPUTATION  INTERVAL  (HOURS)  _  -  2-  _  _ 

-A-JBST-  OA-Y-  C1F  —StMuIA.  C m»L-gEJUilQ - 274  t  t  OCT  701— 

LAST  day  of  simulation  RERinO  304  (31  uCT  70) 

NUMBER  OF  days  in  simulation  period  31  .  . ..  . 


OBSERVATIONS  at  RI 
first  day  flF  STunr 
-LAS  T.  PAY-  -OF  SLUOY^ 
NUMBER  OF  DAYS  In 

v£R  mu  f 

PFR Inn 

.  -  . 

tl.50 

274  (  l 

_ 30-4 _ Ui 

UCT  70) 
lict  /m 

- * - - 

study  otMlf'O 

31 

on.  nF 

rater  duality  parameters  at  river  mile 
number  OF  SIMULATION  points 

11.50 

373 

E  A  X  A  A  AAA  —  .  -  _ 

-  simulation 

VALUES 

(SIMULATED-OBS, ) 

OBSERVED 

parameter 

minimum 

“a  X j mum 

MEAN 

SID. DEV. 

mean  std.dev. 

.  VALUES 

FLOw 

8.9 

37.1 

12.7 

5.9 

TEMP 

9.2 

24.9 

17.7 

3,4 

0,0  0,0 

0 

_ _ AX  Y 

_ 4.0 

U.2  .  . 

_ 9.5- 

.6- 

0,0  0,0 

0 

NH3 

,0  80 

.475 

,3Rl 

.084 

0.000  0,000 

0 

Nl)  J 

.220 

1,215 

1.063 

.166 

O.OOU  0,000 

0 

POa 

.100 

.566 

,440 

.103 

0,000  0.000 

0 

ALKA 

2b.o 

3)  .2 

30.0 

.6 

0.0  0,0 

0 

LOG  COL  I 

2.78 

4.44 

3.03 

.“2 

0.00  0.00 

0 

— TOS  - 

- LQO-, 

Lai.- 

. — 1  r  — 

0,  0, 

A 

PH 

7.3 

7.7 

7.6 

.  i 

o.o  0,6 

0 

ROD 

•  5 

n.t 

2.9 

.7 

0.0  _  0,0 

0 

Parameter 


*ATER  UOalITy  POINTS  EXCEEDING  STaNPARD 

_ STAIDARD  .  _ -NUMBER-  -  -  PERCENT. _ 


TEMP 

OXY 

LOG  CULX 
PM 
PH 


32,2  MAX. 

4,0  MIN, 
3,6'*  “AX, 
6,'l  “IN. 
- 8,5  MAX. 


0 

0,00 

0 

0.00 

52 

13.94 

0 

0.00 

0 

0.00 

TABLE  V-13 

GEORGIA  WATER  QUALITY  STANDARDS 
AND 

OCONEE  RIVER  STREAM  CLASSIFICATION 


DRINKING  WATER  STANDARDS 
(waters  requiring  treatment) 


Col iform 

maximum 

<«nQn  MPN/inn  nl 

Dissolved  Oxygen 

minimum 

4  mg/1 

(warm  water  fish) 

pH 

minimum 

6.0 

maximum 

A.  5 

Temperature 

maximum 

90°F 

NO^  as  Nitrogen 

maximum 

10  mg/1 

FISHING  WATER  STANDARDS 

Col iform 

maximum 

4QOH  MPN/lon  ml 

Dissolved  Oxygen 

minimum 

4  ma/1 

(warm  water  fish) 

pH 

minimum 

6.0 

maximum 

8.5 

Temperature 

maximum 

00°  F 

STREAM  CLASSIFICATION 

North  Oconee 

River 

R.M. 

R.M. 

34.8-23.7 

23.7-12.3 

Drinking  Water 
Fishing  Water 

Middle  Oconee 

River 

R.M. 

31.5-21.8 

Drinking  Water 

Middle  A  Main 

Oconee  River 

R.M. 

21.8-  0.0 

Fishing  Water 

NORTH  OCONEE 


RIVER  MILE 


NORTH  OCONEE 


Figure  V-2.  Existing  Land  Use 
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RIVER  MILE 


Figure  V-4.  Existing  Land  Use 
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Figure  V-5.  Existing  Lond  Use 
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Figure  V-6.  Existing  Land  Use 


TABLE  V-14 

POINT  SOURCE  IMPACTS  ON  THE  NORTH  OCONEE  RIVER 
EXISTING  LAND  USE 


RIVER 

MILE 

PROBABLE 

POLLUTANT 

SOURCE 

PARAMETER 

IMPACT 

- . 

REMARKS 

26.0-23.7 

Unknown 

DO 

•B  mg/1 

minor 

Subbasin  1A  *  2 

Col i form 

9000  MPN/100  ml 

significant 

23.2-21.7 

Subbasin  3 

bod5 

2.5  mg/1 

mi  nor 

nh3 

.2  mg/1 

significant 

no3 

.5 

significant 

PD4 

.2  mg/1 

significant 

19.9-19.3 

STP 

DO 

.6  mq/1 

minor 

BOD5 

4  mg/1 

minor 

nh3 

.6  mg/1 

sianificant 

N03 

1.1  mg/1 

significant 

P04 

.7  mg/1 

significant 

TABLE  V-15 

POINT  SOURCE  IMPACTS  ON  THE  MIODLE  AND  MAIN  OCONEE  RIVER 

EXISTING  LAND  USE 


m 

PROBABLE 

POLLUTANT 

SOURCE 

PARAMETER 

IMPACT 

REMARKS 

22.5-17.0 

STP  *  Sub- 

NH, 

.1  mq/1 

significant 

basin  6A  R  6B 

no3 

.2  mq/1 

significant 

P04 

.1  mn/1 

significant 

BOD5 

1.6  mq/1 

minor 

12.3 

North  Oconee 

NH 

.1  mg/1 

significant 

"4 

.1  mg/1 

significant 

7.2 

Subbasin  16 

BOD5 

.5  mo/1 

minor 

nh3 

.1  mg/1 

significant 
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Analysis  of  Alternative  Future  C 


The  STORM  results  derived  from  alternative  future  C  land  use  condition 
were  imposed  on  the  same  initial  river  quality  condition  as  defined  for 
existing  land  use  in  Table  V-10.  The  tributary  base  flow  quality  condition 
used  for  alternative  C  land  use  is  shown  in  Table  V-16.  Table  V-17  shows 
an  example  of  a  statistical  summary  of  the  water  quality  condition  at  a 
random  point  along  the  river  for  alternative  future  C  land  use  condition. 
The  critical  values  for  some  of  these  parameters  at  various  locations 
along  the  river  have  been  plotted  in  Figures  V-7  to  V-12  to  show  a  river 
water  quality  profile  for  the  most  critical  condition  occurring  during 
the  period  simulated  for  dissolved  oxygen  (DO),  ammonia  (NH3),  nitrate 
(NO-j),  phosphate  (PO^),  log  coliform  bacteria  and  5-day  carbonaceous 
biochemical  oxygen  demand  ( BOD^ ) . 

The  water  quality  condition  for  alternative  future  C  land  use  condition 
exceeds  the  Georgia  State  Water  Quality  Standards  (i.e..  Table  V-13) 
similar  to  the  existing  condition. 

Sample  graphical  results  of  the  simulations  are  shown  in  Appendix  R. 
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TABLE  V-16 

TRIBUTARY  BASE  FLOW  QUALITY  DATA]/ 
ALTERNATIVE  C  LAND  USE 
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TABLE  ¥-17 

WATER  QUALITY  AT  RIVER  MILE  11.5 
ALTERNATIVE  FUTURE  C  LAND  USE 


*»*  ••«**«*««  *»»  *  ******•«.«»*»,*,**  A  ***  AA _ _ _ _ 
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Figure  V-7.  Comparison  of  Water  Quality  Due  to  Changing  Land  Use 
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V-8.  Comparison  of  Water  Quality  Oue  to  Changing  Land  Use 
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Figure  V-9.  Comparison  of  Water  Quality  Due  to  Changing  Land  Use 
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Figure  V-10.  Comparison  of  Water  Quality  Due  to  Changing  Land  Use 


Comparison  of  Water  Quality  Due  to  Changing  Land  Use 
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-12.  Comparison  of  Water  Quality  Due  to  Changing  Land  Use 


Impact  of  Alternative  Future  C 

The  water  quality  Impact  of  alternative  future  C  land  use  is  shown  in 
Figures  V-7  to  V-12.  They  have  been  summarized  in  Table  V-18. 


TABLE  V- 1 8 

WATER  QUALITY  IMPACTS  DUE  TO  CHANGING  FROM 
EXISTING  LAND  USE  TO  ALTERNATIVE  FUTURE  C 


"  MaGnItud 
(mg/1) 

E  LOCATION 

(river  mile) 

S 

BOD, 

1-2 

North  Oconee 

25-12.3 

minor 

D 

1 

Middle  Oconee 

25-0 

minor 

DO 

.2-. 4 

North  Oconee 

22-12.3 

minor 

.2 

Middle  Oconee 

18-15 

minor 

NH, 

.  1-.2 

North  Oconee 

25-12.3 

significant 

.1-.2 

Middle  Oconee 

25-0 

significant 

NO, 

.  1-  .4 

North  Oconee 

25-12.3 

significant 

0 

.1-.3 

Middle  Oconee 

25-0 

significant 

PO, 

.1-.3 

North  Oconee 

25-12.3 

sionificant 

.1-  .2 

Middle  Oconee 

25-0 

significant 

Coliform 

(MPN/lOOmT) 

3,000- 

10,000 

North  Oconee 

33-12.3 

significant 

3,000- 

10,000 

Middle  Oconee 

25-n 

significant 

Remarks  in  Table  V- 1 8  concerning  nutrients  having  sinnificant  impacts  refer 
to  the  potential  impact  on  algal  production  in  non-turbid  water.  Unless 
significant  improvement  occurs  in  the  turbidity  of  the  Oconee  River,  this 
potential  will  not  be  realized. 

In  general,  the  sources  of  the  increased  pollutants  due  to  changing  land  use 
are  the  same  as  those  defined  for  existing  conditions  in  Tables  V-14  and  V-15. 
Concentrations  of  pesticides,  heaw  metals  and  other  parameters  not  speci¬ 
fically  mentioned  were  not  evaluated  in  this  study. 


GRID  CELL  SEDIMENT  TRANSPORT  INVESTIGATIONS 


Introduction 

A  distributed  parameter,  structure  imitating  model  was  developed  for  the 
calculation  of  land  surface  erosion  and  deposition.  The  phenomena 
simulated  in  the  model  are:  rai nfall -runoff ,  runoff  accumulation  and 
distribution,  detachment  of  soil  by  rainfall,  transport  of  detached  soil 
by  runoff,  scour  by  runoff,  and  deposition.  Application  of  the  model  to 
laboratory  test  data  yielded  encouraging  results.  The  application  to  a 
watershed  in  this  study  was  unsuccessful,  however,  because  of  unsatis¬ 
factory  topographical  information. 

Model  Description 

The  model  performs  calculations  on  a  cell-by-cell  basis.  The  direction 
and  velocity  of  runoff  are  determined  from  topographical  information 
Imbedded  In  the  data  base.  A  steady  state  process  Is  assumed.  Details 
of  the  computations  are  given  below. 

Rainfall-Runoff 

The  very  simple  "rational  formula"  was  used  to  generate  runoff  from  each 
cell.  The  runoff  from  any  cell  Is: 

Q  *  ciA 
where : 

0  *  discharge,  cfs 
c  *  runoff  coefficient 
1  *  rainfall  intensity,  in/hr 
A  *  cell  area,  acres 
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The  value  of  the  coefficient,  c,  was  determined  from  a  combination  or 
hydrologic  soil  type  and  land  use  as  shown  in  Table  V-19.  The  values 
in  the  table  are  judgemental  and  have  not  been  calibrated. 

It  is  believed  that  use  of  this  rainfall -runoff  relationship  is  justified 
because  of  the  small  scale  of  the  cells  and  the  steady  state  nature  of 
the  simulation.  The  method  for  accumulating  runoff  from  individual  cells 
Is  described  below. 

Runoff  Accumulation 

Consider  a  typical  cell  (I,  J)  and  its  eight  neighboring  cells  as  shown 
in  Fig.  V-13.  If  any  of  the  neighboring  cells  are  at  higher  elevations, 
a  portion  of  runoff  generated  at  those  cells  will  reach  cell  (I,  J).  The 
runoff  generated  within  cell  (I,  J)  is  found  by  the  rational  formula  and 
added  to  the  sum  of  all  the  contributions  from  higher  cells  to  give  the 
total  discharge  passing  out  of  cell  (I,  J).  This  discharge  Is  evenly 
distributed  among  all  neighboring  cells  of  lower  elevation.  For  this 
reason  computations  must  proceed  from  higher  elevations  to  lower  which 
requires  that  the  data  bank  first  be  sorted  by  elevation.  Note  that, 
since  the  sediment  moves  with  the  runoff,  this  portion  of  the  calculations 
also  determines  the  paths  that  the  sediment  takes.  No  runoff  (or  sediment) 

Is  passed  between  cells  of  equal  elevation. 

Soil  Detachment  by  Rainfall 

Following  a  suggestion  by  Foster  and  Meyer  [14]  it  is  assumed  that  rate  of 
detachment  Is  proportional  to  the  square  of  the  rainfall  Intensity  with 
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the  proportionality  constant  being  the  erodlblllty  factor  (K)  In  the 
universal  soils  loss  equation.  The  formula  used  Is: 

RE  =  a  (K)(Ex)(A)(i2) 
where: 

RE  *  rate  of  soil  detachment  by  rainfall,  tons/hr 

a  *  an  empirical  constant,  the  value  used  (0.0002)  was  based  on 
very  sparse  data 

K  *  soil  erodlblllty  factor  in  the  universal  soils  loss  equation, 
related  to  soil  type 

Ex  =  percent  of  the  cell  area  exposed  to  rainfall,  related  to  land 
use  as  shown  on  Table  V-19 

? 

A  =  cell  area,  ft 
i  *  rainfall  intensity,  in/hr 

In  this  study,  the  rainfall  Is  assumed  to  be  uniformly  distributed  over 
the  basin  and  all  cells  are  assumed  to  be  of  equal  size.  Therefore, 
only  the  values  of  K  and  Ex  change  from  cell  to  cell  to  reflect  spatial 
land  use  variation. 

Sediment  Transport  by  Runoff 

The  hydraulics  of  the  flow  must  be  further  defined  before  sediment  transport 
calculations  can  begin.  To  define  the  hydraulics  certain  Important 
assumptions  must  be  made.  The  formation  of  rills  and  gullies  Is  Important 
to  the  runoff  hydraulics,  however  no  means  of  predicting  their  formation 
or  ultimate  size  was  found.  Therefore,  It  was  assumed  the  runoff  between 
cells  occurs  as  sheet  flow.  The  width  scale,  B,  of  the  runoff  Is,  therefore, 
related  to  the  total  length  of  a  cell  boundary  and  the  number  of  neighboring 
cells  to  which  runoff  Is  passed  (those  with  lower  elevations).  If  flow  goes 
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to  all  eight  neighbors,  B  equals  the  total  length  of  the  cell  coundary; 
if  flow  goes  to  only  one  cell,  B  is  one-eighth  of  the  total,  etc.  From 
the  width,  B,  a  typical  depth  of  flow,  y  for  each  outflow  path  is  calculated 
from  the  Manning  eq.: 

y  ,  Qn  _ 

1.486  B  S  1/2 

where: 

y  »  depth  of  flow,  ft. 

Q  *  total  discharge  from  cell,  cfs 
n  =  Manning's  n 
B  =  width  of  outflow  path,  ft. 

S  *  slope  to  particular  cell,  difference  in  elevation  divided  by 
distance  between  cell  centroids 

Sediment  transport  rates  and  erosion  or  deposition  by  the  runoff  are 
based  on  a  simple  DuBoys  relationship: 

TC  *  CS(t  -  xm)  b 
where : 

TC  *  transport  capacity  of  any  single  outflow  path,  tons/hr 

CS  *  transport  coefficient,  related  to  representative  grain  size: 

CS  3  52.3  where  0  is  grain  size  in  millimeters, 

empirically  related  to  soil  type 

2 

t  3  shear  stress,  Ib/ft 

t  *  yyS,  where  y  3  unit  weight  of  water  and  y  and  S  are  as  previously 
defined 

2 

rm  =  critical  shear  stress  below  which  no  transport  occurs,  Ib/ft 
related  to  grain  size 

b  3  width  of  an  Individual  outflow  path.  Constant  at  one-eighth 
the  total  length  of  a  cell  boundary 
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Whether  scour  or  deposition  occurs  alonq  any  particular  outflow  path  depends 
upon  whether  the  sediment  load  is  less  than  or  greater  than  the  transport 
capacity.  The  load  to  any  outflow  nath  is  calculated  as  follows: 
g  *  (GI  +  RE)  (Q0PS/O0) 
where: 

g  *  sediment  inflow  to  any  particular  outflow  path,  tons/hr 

GI  =  total  rate  of  sediment  inflow  to  the  cell  from  neiqhborinq  cells 
of  higher  elevation,  tons/hr 

RE  3  rate  of  sediment  detachment  by  raindrop  within  the  cell,  ton/hr. 

Q0PS  3  runoff  following  individual  outflow  path,  cfs,  equal  to  Q0  divided 
by  the  number  of  outflow  paths 

QP  3  total  runoff  from  cell,  cfs 

If  g  is  less  than  the  transport  capacity,  erosion  occurs  along  the  outflow 
path.  The  actual  transport  rate  for  that  path  is  calculated  by  the 
following: 

G2N  3  q  (1  -  Ex)  +  ((1  -  F)  q  +  F(TC)Ex) 

where: 

G2N  3  sediment  transported  out  of  cell  alonq  anv  given  outflow  path, 
tons/hr 

F  =  a  relaxation  factor  if  F  =  0,  outflow  =  inflow;  if  F  =  1, 
outflow  3  transport  capacity 

Other  symbols  are  as  previously  defined.  The  exposure  factor  appears 
because  erosion  can  only  occur  where  the  soil  is  available. 


If  g  is  greater  than  the  transport  capacity,  deposition  occurs  and  the 
outflowing  load  is  calculated  as  follows: 


G2N  «  (1  -  F)  g  +  F  (TC) 


All  symbols  have  been  defined.  The  exposure  does  not  appear  because 
deposition  can  occur  everywhere. 

These  values  are  added  to  the  inflowing  load  (Gl)  of  the  neighboring 
cells. 
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Application  to  Laboratory  Data 

The  algorithm  for  calculation  of  land  surface  erosion  and  deposition 
described  above  was  tested  by  comparing  calculated  erosion  rates  with 
those  measured  In  a  laboratory.  While  such  a  test  does  not  constitute 
rigorous  verification,  it  can  be  used  to  evaluate  the  general  validity 
of  the  approach  and  Identify  some  inadequacies. 

The  laboratory  test  data  used  [15]  Is  from  a  5  foot  by  16  foot  plot. 

It  was  modeled  using  25-1  by  1.6  feet  rectangular  cells.  The  rough¬ 
ness  value  (n  a  0.022)  and  runoff  coefficient  (C  *  0.97)  were  based  on 
measurements  made  during  the  experiment.  Although  several  slopes  were 
tested,  only  the  10%  slope  condition  was  modeled.  This  slope  was  re¬ 
flected  In  elevations  assigned  to  the  Individual  cells.  The  simulated 
rainfall  Intensities  were  used  in  the  program  and  the  calculated  weight 
of  sediment  transported  to  the  bottom  of  the  plot  compared  with  that 
measured.  The  results  are  shown  in  Table  V-20. 

The  same  set  of  coefficients  was  used  for  all  rainfall  intensities  In¬ 
dicating  that  the  functional  relationship  between  rainfall  and  erosion 
used  Is  reasonable.  A  relaxation  coefficient,  F,  of  zero  had  to  be  used. 

A  zero  F  indicates  that  the  runoff  has  sufficient  transport  capacity  to 
carry  all  the  sediment  produced  by  raindrop  erosion.  Apparently  the 
transport  capacity  calculated  by  the  DuBoys  relation  was  too  small. 

This  could  be  due  to  the  assumed  runoff  hydraulics  being  Inadequate, 
or  the  critical  shear  stress  not  being  appropriate  for  land  surface 
erosion.  The  critical  shear  stress  was  taken  from  a  Shield's  diagram  [16] 
developed  for  open  channel  flow. 
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TABLE  V-20 


Rainfall 

Intensity 

(In/hrj 

Measured 

Wt.  of  Sediment 
(lbs/hr) 

Calculated 

Vft.  of  Sediment 
(Ibs/hr) 

1.25 

5.3 

7.6 

2.25 

27.1 

24  i 

i 

3.65 

67 

66 

4.60 
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Application  to  Sandy  Creek  Basin 

The  model  was  applied  to  the  lower  16  square  miles  of  the  Sandy  Creek 
watershed.  An  existing  detailed  data  bank  was  available  for  this  area. 

The  spatial  variation  of  various  parameters  within  the  watershed  was 
described  using  a  total  of  9208  grid  cells.  The  variables  used  from  the 
data  bank  were:  cell  elevation,  soil  type,  hydrologic  soil  type  and  land 
use.  A  topographic  map  of  the  basin  is  shown  in  Figure  V-14. 

Runoff  coefficients  were  determined  from  land  use  and  hydrologic  soil 
type  as  given  in  Table  V-19.  Also  shown  are  the  exposure  values  assigned 
to  the  various  land  uses.  Descriptions  of  the  land  uses  are  given  in 
Chapter  III,  and  soil  erodlbility  factors  (K)  in  Table  V-5.  After  several 
runs  and  mapping  of  computed  discharges  and  sediment  loads,  a  basic  data 
problem  was  identified  which  prevented  completion  of  the  apolicatlon. 

This  problem  Is  discussed  in  detail  below. 

Topographic  Data  Problems 

The  procedure  used  relies  on  topographic  data  in  the  form  of  an  elevation 
for  each  cell.  Differential  elevations  between  neighboring  cells  drive  the 
runoff  calculation.  If  the  cell  elevations  are  truly  representative  of  the 
topography  of  a  drainage  basin,  every  cell  but  one  will  have  at  least  one 
outflow  path.  The  exception  is  the  outlet  on  the  watershed  boundary. 

The  cell  elevations  In  the  Oconee  study  were  manually  assigned  from  a  base 
topographic  map.  This  procedure  resulted  in  many  cells  for  which  no  outflow 
path  existed.  Out  of  9208  total  cells,  334  had  no  outflow  path.  Since  no 
runoff  or  sediment  can  pass  through  a  cell  with  no  outflow  path,  the 
calculations  for  all  downhill  cells  are  erroneous. 
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Several  "smoothing"  algorithms  were  tried  to  insure  that  all  cells  had 
outflow  paths.  Of  these,  the  best  appeared  to  be  the  following:  set  the 
no  outflow  cell's  elevation  equal  to  the  average  elevation  of  the  lowest 
and  second  lowest  neighboring  cells.  This  guarantees  that  the  cell  in 
question  has  an  outflow  path,  but  may  eliminate  the  only  outflow  path  a 
neighboring  cell  has.  The  algorithm  worked  fairly  well;  after  several 
passes  the  number  of  no  outflow  cells  was  reduced  from  334  to  about  18, 
which  could  not  be  further  reduced  by  successive  application  of  the 
algorithm.  Unfortunately,  this  was  still  too  high  to  yield  acceptable 
results,  as  shown  on  Figure  V-15.  Happed  on  that  Figure  are  cell-by-cell 
discharges  generated  by  1-1nch  of  rain  in  one  hour.  The  tendency  for  the 
runoff  to  accumulate  In  stream  channels  can  readily  be  seen.  Note  also, 
however,  the  effect  of  cells  with  no  outflow  path  (indicated  by  large 
dots).  The  effect  Is  to  disconnect  the  basin,  so  that  runoff  from  upper 
portions  of  the  basin  does  not  pass  through.  Sediment  transport  calculations 
based  on  this  runoff  pattern  are,  of  course,  meaningless. 

If  a  workable  algorithm  cannot  be  developed  for  editing  the  topographic 
data,  the  procedure  for  initially  determining  cells'  elevations  must  be 
modified.  Automatic  interpolation  is  one  possibility.  This  process 
would  rarely  produce  neighboring  cells  of  exactly  the  same  elevation. 

Another  procedure  being  investigated  describes  topography  as  an  array 
of  triangular  elements.  Elevations  are  prescribed  at  the  vertices  and 
vary  linearly  within  each  element. 
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APPENDIX  A 

SMITH  REPORT 


ADAPTION  OF  WATER  QUALITY  -  ECOLOGICAL  MODEL 


TO  THE  OCONEE  RIVER  SYSTEM 


By 

Donald  J.  Smith 
Tetra  Tech,  Inc. 

3700  Mt.  Diablo  Boulevard 
Lafayette,  California  94549 


BACKGROUND  AND  PURPOSE 

The  U.S.  Army  Corps  of  Engineers  Hydrologic  Engineering  Center 
(HEC)  is  adapting  a  dynamic  water  quality-ecological  model  to  a  portion 
of  the  Oconee  River  in  Georgia.  While  the  model  is  designed  to  calcu¬ 
late  the  population  dynamics  of  algae,  zooplankton,  benthic  animals  and 
fish,  detailed  calibration  of  that  portion  of  the  model  is  beyond  the 
limited  scope  of  this  project.  Since  the  primary  purpose  of  the  model 
will  be  to  evaluate  the  transient  water  quality  impact  of  storm  runoff 
and  waste  discharges,  the  organic  sediment  and  the  biological  section 
of  the  model  will  remain  constant  during  the  simulation. 

In  lieu  of  modeling  those  parameters  which  were  held  constant, 
pertinent  reports  of  water  quality  and  biological  surveys  were  reviewed 
to  estimate  their  values. 

This  brief  report  documents  the  findings  of  this  review. 


RIVER  SYSTEM 

The  study  is  limited  to  the  upper  reaches  of  the  Oconee  River  sys¬ 
tem  near  Athens.  Included  is  the  Oconee  River  between  the  Barnett 
Shoals  Dam  and  the  confluence  of  the  Middle  Oconee  and  North  Oconee 
Rivers  and  approximately  twenty- five  (25)  miles  of  the  Middle  Oconee 
River  and  twenty  (20)  miles  of  the  North  Oconee  River. 
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The  Middle  and  North  Oconee  Rivers  are  typically  1/2  to  5  feet 
deep  and  from  50  to  100  feet  wide.  The  Oconee  is  also  typically  50  to 
100  feet  wide  with  depths  up  to  8  feet.  The  average  gradient  is  approxi¬ 
mately  4  feet  per  mile  and  velocities  are  characterized  as  slow  to 
modera te. 

Urban  development  is  limited  to  the  Athens  area.  The  remaining 
watershed  is  rural  with  many  wooded  areas.  A  total  of  approximately 
seven  (7)  million  gallons  a  day  of  municipal  and  industrial  waste  water 
is  discharged  to  the  Middle  Oconee  and  North  Oconee  after  secondary 
treatment.  During  periods  of  high  runoff,  significant  amounts  of  or¬ 
ganic  detritus  and  sediment  are  washed  into  these  rivers  from  the 
watershed. 


WATER  QUALITY 

For  purposes  of  characterizing  the  water  quality,  the  river  system 
can  be  divided  into  two  sections.  Section  One  includes  the  Oconee  River 
and  those  portions  of  the  Middle  Oconee  and  North  Oconee  below  the  two 
Athens  sewage  treatment  plant  outfalls.  Water  quality  in  this  section 
is  influenced  by  the  Athens  sewage  treatment  plant  effluent.  Dissolved 
oxygen  is  lower  and  plant  nutrients,  BOD,  and  total  organic  carbon  are 
higher  than  in  Section  Two.  Water  quality  in  Section  Two,  the  remaining 
portion  of  the  study  area,  is  reasonably  good.  The  water  of  both  sec¬ 
tions  is  quite  turbid  during  periods  of  high  flow.  Levels  of  selected 
water  quality  parameters  reported  by  state  and  federal  agencies  are 
summarized  in  Table  1. 


ALGAE  AND  ZOOPLANKTON 

No  chlorophyll  a  data  or  other  direct  measurement  of  suspended  al¬ 
gae  are  available.  Some  attached  algae  (periphyton)  and  the  macrophyte 
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Podostemum  have  been  observed  where  suitable  rock  substrate  is  avail- 

p 

able.  Nelson  (1962)  reports  Podostemum  levels  of  10  to  15  g/m  (dry 

weight)  on  the  Middle  Oconee  where  ideal  substrate  conditions  exist  and 

none  on  sand  or  mud  substrates.  Suitable  substrate  (bedrock  and  cobbles) 

should  exist  where  velocities  are  sufficiently  high  to  prevent  deposition 

of  sand  and  silt.  If  we  assume  velocities  are  sufficiently  high  with 

bottom  slopes  of  1  foot  in  200  feet  and  average  bottom  slope  of  4  feet 

in  one  mile,  approximately  15%  of  the  substrate  is  suitable.  Fifteen 

percent  of  the  densities  reported  by  Nelson  yield  average  macrophyte 

2 

and  periphyton  densities  from  1.5  to  2.25  g/m  . 

An  examination  of  dissolved  oxygen  data  indicates  that  algal 
photosynthetic  oxygen  production  is  not  significant.  Dissolved  oxygen 
never  exceeded  saturation  and  no  diurnal  variation  was  evident.  No  de¬ 
pletion  of  plant  nutrients  was  observed  in  the  data.  Nelson  also  re¬ 
ported  that  no  detectable  differences  in  dissolved  oxygen  was  observed 
between  upstream  and  below  his  study  area. 

For  modeling  purposes,  both  algae  and  zooplankton  concentrations 
can  be  presumed  low,  near  zero. 


BENTHIC  ANIMALS 

The  make  up  of  the  benthic  animal  (micro invertebrates)  population 
has  been  studied  by  submerging  limestone  substrate  (LSS)  in  the  water 
for  two  months.  The  results  of  these  studies  indicate  that  benthic 
animals  can  survive  throughout  the  study  area  if  suitable  substrates 
exist.  Unfortunately  the  test  results  do  not  include  total  biomass  or 
LSS  surface  areas,  therefore,  population  densities  cannot  be  determined. 

Nelson  (1962)  observed  benthic  animal  population  densities  of  2  to 

10  g/m  (dry  weight)  with  ideal  substrate  conditions.  If  we  also  assume 

fifteen  percent  (15%)  of  the  natural  substrate  is  suitable  for  benthic 

2 

animals,  average  densities  of  .3  to  1.5  g/in  can  be  expected. 
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FISH 


Fish  population  data  within  the  study  area  are  limited  to  one  sam¬ 
pling  event  in  June  1959.  Total  fish  mass  was  reported.  However,  the 
length  of  stream  sampled  was  omitted  making  it  impossible  to  calculate 
fish  biomass  per  mile.  All  fish  collected  were  warm  water  species. 
Approximately  fifty  percent  ( 50%)  were  bottom  feeding  fish. 

Streams  of  this  type  typically  have  a  fish  population  of  100  to 
300  lb/acre  (wet  weight).  Assuming  an  average  channel  width  of 
75  feet,  approximately  600  feet  of  channel  has  a  surface  area  of  one 
acre.  Converting  to  dry  weight  per  mile,  a  total  biomass  of  90  to  260 
lbs/mile  or  40  to  120  kg/mile  is  obtained. 

DETRITUS  AND  ORGANIC  SEDIMENT 

Total  organic  carbon  in  the  water  ranged  from  2  to  7  mg/1  C.  The 
detritus  concentration  in  the  water  is  generally  twice  the  organic 
carbon  level  or  4  to  14  mg/1 .  The  detritus  level  is  generally  a  function 
of  flow  rate,  increased  detritus  occuring  with  increased  flow.  Nelson 
(1962)  attempted  to  correlate  river  discharge  with  detritus  volatile 
solids)  with  some  success.  He  typically  measured  volatile  solids  of 
1  to  10  mg/1  at  moderate  flows  and  up  to  50  mg/1  during  high  flow  periods. 

Nelson  also  measured  settleable  plant  and  animal  detritus  and 
reported  typical  values  of  12  to  20  g/m^.  Below  the  Athens  STP  outfalls, 
suspended  detritus  and  organic  sediment  can  be  expected  to  increase.  A 
10  to  20  percent  increase  in  the  above  values  seems  appropriate. 


Nelson,  Daniel  J.,  and  Scott,  Donald  C.,  1962,  Role  of  detritus  in  the 
productivity  of  a  rock-outcrop  community  in  a  piedmont  stream.  Depart¬ 
ment  of  Zoology,  University  of  Georgia. 
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